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On  the  Generation  of  Stress  end  Deformation 
in  Elastic  Solids  by  High  Powered  Lasers 

Peter  J.  Torvik 


1 .  INTRODUCTION 


The  absorption  of  some  fraction  of  the  energy  in  a  high  powered  laser 
beam  at  the  surface  of  a  solid  generates  stress  and  deformation  through 
several  mechanisms.  These  may  be  broadly  divided  into  two  categories: 

(1)  where  the  stress  arises  through  the  thermal  expansion  of  the  target 
material,  and  (2)  where  the  stress  in  the  target  arises  through  a  pressure 
loading  resulting  from  the  thermal  expansion  of  the  medium  (typically  air). 

The  blow-off  phenomenon,  in  which  the  stress  is  generated  in  the  target  by 
the  rapid  expansion  of  vaporised  target  material,  may  be  regarded  as  a  third 
category  or  may  be  arbitrarily  assigned  to  either  of  the  first  two. 

A  second  rationale  for  classifying  the  various  regimes  of  stress  gener¬ 
ation  is  furnished  by  the  time  scales  of  the  problem.  First,  the  time  scale 
introduced  through  mechanical  inertia  of  the  target  must  be  considered.  If 
the  heating  rates  are  sufficiently  low,  it  is  evident  that  acceleration  terms 
in  the  target  response  can  be  neglected  and  a  static  stress  analysis  per¬ 
formed.  Such  problems  may  be  termed  thermal  stress  problems.  On  the  other 
hand,  under  very  short  heating  times,  significant  temperature  changes  may 
occur  before  significant  deformation  has  occurred.  In  such  cases,  the  influ¬ 
ence  of  target  material  inertia  is  to  serve  as  a  constraint  thereby  rairing 
the  levels  of  thermal  stress.  These  may  be  labeled  as  thermal  shock  problems. 

The  finite  values  of  target  conductivity  introduce  a  second  time  scale 
into  the  target  response.  If  heating  rates  are  very  low,  a  state  of  thermal 
equilibrium  will  result  with  input  energy  being  balanced  by  "losses"  to  the 
surrounding  media  or  through  radial  conductions  away  from  the  heated  area.  1 
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These  may  be  regarded  as  static  thermal  stress  problems .  For  more  rapid 
heating  rates,  changes  in  temperature  must  be  taken  Into  account,  but  the 
stress  and  deformation  analysis  may  be  performed  by  neglecting  mechanical 
inertia  and  performing  a  static  stress  and  deformation  analysis  with  time 
regarded  as  a  parameter.  Such  problems,  where  thermal  inertia  is  considered 
but  mechanical  inertia  is  not,  will  be  termed  quasi-static  thermal  stress 
problems.  In  comparison,  the  thermal  shock  regime  is  seen  to  exist  when  both 
mechanical  inertia  and  thermal  inertia  must  be  considered.  The  fourth  possi¬ 
bility,  that  thermal  inertia  may  be  neglected  but  mechanical  inertia  must  be 
considered,  occurs  in  a  partially  transparent  poor  conductor  where  absorption 
of  a  pulse  occurs  in  depth.  The  initial  temperature  profile  and,  hence,  the 
gradients  and  resulting  stresses  will  form  and  be  propagated  quite  indepen¬ 
dently  of  the  thermal  diffusion  process.  For  purposes  of  classification, 
these  may  be  called  isothermal  shock  problems  and  are  of  particular  signifi¬ 
cance  in  the  afcs<*rgtion  of  x  rays.  However,  this  class  would  appear  to  be  of 
an  academic  interest  only  in  targets  essentially  opaque  to  laser  radiation, 
where  absorption  occurs  at  the  surface  and  thermal  stresses  are  generated 
after  a  temperature  change  arises  through  conduction. 

A  further  classification  of  thermal  stress  problems  may  be  made  into 
problems  of  coupled  and  uncouoled  thermoelasticity.  In  coupled  thermoelas¬ 
ticity,  the  temperature  change  resulting  from  a  mechanical  strain  induced 
adiabatically  is  taken  into  account.  This  effect  is  not  to  be  expected  to  be 
of  significance  in  the  presence  of  the  very  large  temperature  changes  arising 
from  the  absorption  of  intense  laser  radiation.  A  second  deviation  from  clas¬ 
sical  thermoelasticity  arises  if  non-Fourier  solids  are  considered  in  which 
the  heat  conduction  law  is  modified  so  as  to  provide  finite,  rather  than  infi¬ 
nite,  velocities  for  the  propagation  of  thermal  disturbances.  Solutions  to 
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such  nonclassical  problems  should  be  reviewed  In  order  to  determine  if  this  ef¬ 
fect  is  of  importance  in  the  generation  of  thermal  stresses  by  laser  heating. 

Three  of  these  four  regimes  of  uncoupled  thermoelasticity  will  be  consid¬ 
ered  in  the  second  section  of  this  monograph  with  the  goals  being  as  follows: 
(1)  to  determine  circumstances  under  which  significant  thermal  stresses  will 
be  generated,  and  (2)  to  provide  readily  used  formulae,  charts  and  algorithms 
for  estimating  the  magnitude  of  the  thermal  stress.  Consideration  of  isother¬ 
mal  shock  will  be  left  for  a  later  work. 

In  the  third  section,  the  generation  of  stress  and  deformation  by  mecha¬ 
nisms  involving  a  pressure  loading  will  be  considered  with  the  emphasis  again 
being  on  determining  circumstances  under  which  the  stresses  and  deformations 
are  significant  and  on  providing  means  of  estimating  the  magnitude  of  the 
effect. 

In  the  fourth  section,  damage  mechanisms  will  be  treated.  Of  particular 
Interest  is  determining  the  likelihood  that  fluctuations  in  stress  by  what¬ 
ever  mechanism  produced  will  be  of  such  magnitude  as  to  cause  failure  through 
catastrophic  (unstable)  crack  growth  or  eventual  fracture  through  crack  prop¬ 
agation.  In  order  to  better  assess  the  significance  of  stresses  and  deforma¬ 
tions  produced  by  laser  heating,  comparisons  will  be  made  with  other  damage 
mechanisms,  notably  melt  through  and  thermal  softening. 


II.  UNCOUPLED  THERMOELASTICITY 

The  classification  scheme  introduced  in  the  Introduction  was  based  on  two 
time  scales.  The  first  was  the  thermal  time  scale,  which  arises  from  the  dif¬ 
fusive  nature  of  the  heat  conduction  equation.  As  will  be  seen  in  more  detail 
later,  solutions  to  one-dimensional,  linear  heat  conduction  problems  typically 
show  a  dependence  on  time  and  spatial  coordinate  in  the  combination  xt/x*, 
where  k  is  the  thermal  diffusivlty,  t  is  time  and  x  the  characteristic 

coordinate.  Thus,  for  times  t  long  compared  to  a  critical  time 

tc  =  L*/k  (1) 

where  L  is  a  characteristic  dimension  of  the  object,  the  temperature  field 
may  be  regarded  as  static.  Conversely,  for  t  <  t^,  significant  changes  in 
the  temperature  field  will  occur  and  must  be  accounted  for  through  solving 
the  diffusion  equation  for  the  temperature  field  subject  to  appropriate 
thermal  boundary  conditions. 

The  second  time  scale  of  interest  is  the  propagation  time  for  elastic 
disturbances 


L/C 


e 


where  L  is  again  a  critical  dimension  of  the  heated  object  and  Cft  is  an  elas¬ 
tic  wave  speed, of  order  2-5  km/sec  in  an  unbounded  solid.  Table  I  depicts 
the  result  of  this  classification  and  sets  down  the  terminology  to  be  used 


in  this  section. 


TABLE  I 


Classification  of  Uncoupled  The rmoe last ic 

Thermal  Inertia 
t  <  L1/*  t  >  L*/k 

Static  Thermal  Stress  X 

Quasi  Static  Thermal  Stress  X 

Isothermal  Shock  X 

Thermal  Shock  X 


Problems 


Mechanical  Inertia 

t  <  L/C  t  >  L/C 
e 

X 

X 


X 

X 


For  an  isotropic  Fourier  solid,  the  components  of  the  heat  flux  vector, 
q^,  are  related  to  the  spatial  gradients  of  the  temperature  field  through 


q 


i 


(2) 


and  conservation  of  energy  leads  to 


p£  = 


VtJ  *  qi,i 


(3) 


where  k  is  the  thermal  conductivity,  p  is  the  density,  6  is  the  internal  ener¬ 
gy,  o .  .  are  the  stresses,  and  c. .  the  strains*  For  an  elastic  solid 


where  the  entropy  n  is  given  by 


♦  =  £(Eij(  T)  -  n(cij,  T)T 


is  the  free  energy.  From  (3),  (4)  and  (6),  Equation  (5)  follows,  and  that 


°U  *  *  ^ 


The  specific  heat  at  constant  volume,  C^,  is  defined  from  Equation  (3),  with 

*u  - 0>  or 


pCvT  -  “£  •  -"i,I 


Equation  (4)  may  be  written  as 


_q  =pTf_3n_;  +  m  i  "1 

qi,i  pi  I  3ei j  eij  +  3T  | 

=  pT  r..-alL.;  .Hi  i] 

pl  j  3CljaTci.5  W  | 


substituting  Equation  (7)  and  noting  from  (8)  and  (10)  that 


we  arrive  at 


pC  =  -pT 
v  3T* 


-qi,i  =  T(-"3T  )eij  +  PV 
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For  an  elastic  material  undergoing  temperature  change 


°ij  "  XSijGkk  +  2weij  “  ^3X  +  “  V  ^13* 

where  X,p  are  the  Lame  constants,  «T  is  the  coefficient  of  thermal  expansion, 

and  T  is  the  uniform  temperature  in  the  reference  (unstrained)  state.  Sub- 
o 

stituting  (13)  and  (2)  into  (12)  leads  to 

k  irk =  pcv  It  +  *<3X  ♦  ^vi*  (u) 

i  't 

Linearization  (by  fixing  T  at  a  constant  value)  leads  to  the  coupled  theory 
of  thermoelasticity.  Neglecting  the  last  term  leads  to  the  familiar  governing 
differential  equation  of  linear  heat  conduction 

KV*T  =  |~  (15) 

where  the  thermal  diffusivity,  k,  is  defined  by 

k  =  —■  (16) 

PCV 

This  development  is  a  summary  of  the  more  comprehensive  treatment  given  in 
standard  sources,  such  as  Boley  and  Weiner^-.  Numerous  solutions  to  Equa¬ 
tion  (15)  for  various  boundary  conditions  are  given  in  the  standard  work  on 

2 

heat  conduction  in  solids  .  Solutions  of  particular  Interest  to  laser  heating 

3 

problems  have  also  been  identified  . 

Boley ,  B.  A.  and  Weiner,  J.  H.,  Theory  of  Thermal  Stresses,  John  Wiley, 
New  York,  1960. 

2 

Carslaw,  H.  S.  and  Jaeger,  J.  C.,  Conduction  of  Heat  in  Solids,  Oxford 
University  Press,  1959. 

3 

Torvik,  P.  J.,  Thermal  Response  Calculations  and  Their  Role  in  the 
Design  of  Experiments,  AFIT  TR  73-6,  1973, 
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The 


strains  are  related  to  displacements,  U^,  through 


1  ,aui  IV 
cij  =  2  <8^  +  TT^ 


(17) 


in  the  infinitesimal  theory.  Conservation  of  linear  momentum  requires  that 


8o 

— ii  =  PU 

3x  p  i 


(18) 


where  the  stresses  o^  are  related  to  strains  through  Equation  (13)  and, 
hence,  to  displacements.  Substitution  of  (13)  and  (17)  into  (18)  produces 


,  aHv 

<A  *  ->  dr  <ir>  *  “ 


3x^8Xj 


-  <3>  +  2p>  “t  TT 


3*Ui 

"  p  Tt7* 


(19) 


The  problem  of  uncoupled  thermoelasticity  is  the  solution  of  Equation  (19) 
where  the  temperature  gradients  are  determined  from  the  solution  of  (15) 
subject  to  appropriate  boundary  conditions. 
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A.  Static  Thermal  Stresses 

Thermal  stress  arising  from  a  static  temperature  field  will  not  be  en¬ 
countered  frequently  with  high  levels  of  Incident  intensity,  but  one  partic¬ 
ular  case  is  of  such  significance  as  to  merit  consideration. 

Consider  a  flat  sheet  of  f.iickness  1  and  large  lateral  dimension.  Energy 
is  absorbed  at  late  I  per  unit  area,  per  unit  time,  uniformly  over  a  circle 
of  radius  a.  We  assume  that  the  temperature  under  the  beam,  T^,  is  essen¬ 
tially  uniform  and  steady.  An  energy  balance  on  the  disk  of  thickness  1  and 
radius  a  immediately  under  the  beam  leads  to 

I  ita1  =  2(*a*)[hT1  +  eo{(T  +  T.)"  -  T  *}]  -  2walk  —  I  (1) 

a  1  o  1  o  J  3r  j 

*  r=a 

The  first  term  on  the  right  represents  losses  to  the  surroundings  at  tempera¬ 
ture  T  ,  and  the  second,  losses  from  the  nested  disk  due  to  radial  conduction, 
o 

For  simplicity,  the  loss  term  may  be  approximated  by 

h*TL  =  (h  +  cof*)  Tl  (2) 

where 

T*  =  [(T  +  T,)1  h  T  *1[2T  +  T, ]  (3) 

o  1  oJ*-olJ 

Tj  being  the  maximum  temperature  reached  relative  Vo  the  surroundings,  e  being 
the  emissivity,  o  being  the  Stefan-Boltzman  constant,  and  T^  being  the  final 
temperature  of  the  heated  disk  measured  with  respect  to  the  temperature  of  the 
surroundings. 

The  approximation  of  a  uniform  temperature  under  the  heated  spot  is  ap¬ 
propriate  from  two  points  of  view.  First,  the  precise  details  of  the  temper- 
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•cure  distribution  are  dependent  on  local  properties  of  the  bean,  generally 
not  well  characterised.  Secondly,  local  surface  neltlng  under  the  bean  will 
serve  to  "smear"  the  tenperature  into  a  uniform  value,  namely,  the  neltlng 
temperature.  It  is  presumed  in  Equation  (1)  that  the  thermal  inertia  is 
negligible,  i.e. 

pC  ta* IT  «  I  ta*  (4) 

P  * 

Outside  the  heated  region,  r  _>  a,  we  have 

kV*T  -  2h*T/l  =  pC  T  (5) 

P 

with  boundary  conditions  T(r)  =  T. ,  T(»)  =  0,  and  the  gradient  at  r  =  a  given 
by  Equation  (1).  Again  neglecting  the  thermal  inertia  and  defining  a  Biot 
nunber 


6  * 


h*l 

k 


(6) 


we  find 


T  =  AK  (/2B  r/i) 
0 


(7) 


where 


A  =  T.(K  (/2B  a/1))'1 
1  o 


(8) 


Here  K  is  the  modified  Bessel  function  of  the  second  kind  and  order  sero.  If 
o 

T.  <  T  ,  the  melting  temperature,  we  see  from  Equation  (1)  that 
1  ""  n 


.  __  K  (/2B  a/1) 

1  <  2h*T  +  2  -  /2B  v  ,7*r — 777  T. 

a  -  1  a  •Z1)  1 


(9) 
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We  a«y  expect  losses  from  the  heated  area  to  the  surroundings  (the  first  tern 
on  the  right  hand  side)  to  be  negligible  compared  to  losses  due  to  radial  con¬ 
duction  (the  second  term)*  For  small  arguments, 


K  (z)  ~  -lnz  (10a) 

o 

Kl(z)  -  1/z  (10b) 


Thus 


T1  ~ 


I  a1 

a 

2k  1 


In  (- 


a/2B 


when  /2fi  y  <<  1 


(ID 


and  the  temperature  for  small  r  decays  logarithmically. 
From  the  temperature  distribution 


r  <  a 

(12a) 

K  (/2B  r/1) 

— -  ,  r  >  a 

(12b) 

K  (/2B  a/1) 
o 

with  given  by  (9)  or  (11),  as  appropriate,  we  uay  compute  thermal  stresses 
using  results  given  previously.^ 

For  any  axially  symmetric  temperature  distribution  T(r),  the  stresses  and 
displacements  are: 


rr 


ee 


,-.EfrTH  EL2 

r*  r  r  +  1-v  ( l+v)r* 

-sf 


EC1  ECj 

Trdr  -  *  (Uv)rT 


(13) 


(14) 


Boley,  B.  A.  and  Weiner,  J.  H.,  Theory  of  Thermal  Stresses,  John  Wiley, 
New  York,  1960,  pp  288-290. 
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u 

r 


utvy» 

r 


Trdr  +  C,  r  +  — 
1  r 


(1 


for  Arbitrary  C.  Here  we  require: 

In  Region  1:  (r  <  a)  that  remain  bounded. 

in  Region  17:  r  >  a  that  a  -  0  as  r  +  00. 

rr 

on  r  =  a  that  the  normal  traction  and  displacements  of  Region  I  match 
those  of  Region  II.  If  we  neglect  the  change  in  modulus  resulting  from  the 
temperature  increase ,  for  r  _<  a: 


o 


rr 


-EaT^/2 

-EoTj^/2 

(1+v)aT1r/2 


( 

( 

( 


for  r  >  a: 


a 

rr 


q( 1+v ) 
r 


rTdr  - 


ET^a* 

2r* 


ET.aa* 

rTdr  +  — jp - EaT 


rTdr  :•  (1+v) 


oT^a 

2r 


( 


( 


( 


For  the  temperature  field  given  by  Equations  (7)  and  (8),  for  r  >  a 


1 


K  fori 


[p^K.Cp*)  -  prK.  (pr)] 


rTdr  = 


Approximate  solutions  then  give  the  following  results 


M 

r*/a 

5 

1.2 

10 

1.08 

20 

1.04 

From  these  few  values,  we  see  that  the  ratio  r*/a  does  indeed  approach 
unity  for  large  pa,  as  expected.  Thus,  the  region  of  large  compressive  cir¬ 
cumferential  stress  becomes  small,  the  radius  at  which  tensile  stresses  de¬ 
velop  is  reduced,  and  the  magnitude  of  those  tensile  stresses  increases  as  pa 
is  increased.  Recalling  that  a  large  pa  corresponds  to  a  poor  conductor,  and 
that  stresses  decay  as  r  increases,  we  see  that  large  tensile  circumferential 
stresses  are  a  peculiar  problem  to  poor  thermal  conductors. 


One  aspect  of  these  results  is  somewhat  astonishing.  We  note  from 
Equations  (16)  and  (17)  that  the  stresses  under  the  heated  region  are  inde¬ 
pendent  of  the  value  of  pa  and  hence  of  the  quantity: 

Jli.il  (28) 

k  1 

l.e,,  the  thermal  conductivity  outside  the  heated  region  does  not  affect 
stresses  within.  Thus,  in  good  conductors,  such  as  ductile  metals,  large 
maximum  shear  stresses  will  be  found  under  the  heated  region,  but  tensile 
stresses  outside  will  be  small.  Thus  failure  will  be  localised  to  the  heated 
region.  For  poor  conductors,  however,  large  tensile  stresses  will  also  be 
found  outside  the  heated  region,  thus  introducing  the  possibility  of  cracking, 
or  crack  growth,  in  the  portions  of  a  target  .iot  heated  directly. 

For  computational  purposes,  let  us  define 

F^(pa)  =  paKj(pa)  +  fop  Kp(pa)  (29) 

Fj(pr)  =  prK^(pr)  +  (pr)*Ko(pr)  (30) 

Then  Equation  (:*6)  may  be  rewritten 


°ae  ■  F2<Pr)} 

EoT,  =  (pa)*K  (piO 
1  o 


Figure  1  depicts  graph!: si ly  the  functions  F^,  for  a  range  of  arguments 
and  shows  tl  at  small  pa  does  not  give  rise  to  small  pr.  Rather,  we  note  that 
the  radius,  r*,  at  which  the  circumferential  stress  becomes  tensile  remains 
large.  From  Equation  (31)  tensile  stresses  develop  for  F^(pa)  =  F^pr*). 
From  the  graph  we  may  deduce  the  following 
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HajKVg"  ■  i  TT  na/i 


£• 

£T* 

r»/s 

.01 

1.2 

120 

.1 

1.2 

12 

.5 

1.4 

7 

1.0 

1.75 

1.75 

2.0 

2.7 

1.35 

5.0 

5.7 

-1.14 

Because  of  the  decay  of  stress  with  r*  (Equation  (31)),  the  tensile  stress 
which  eventually  develops  In  a  good  conductor  will  be  negligible.  Thus  ten¬ 
sile  stresses  of  significance,  i.e.  near  to  the  heated  region,  will  develop 
only  for  poor  conductors.  (It  also  shows  the  asymptotic  expansions  used 
(Equation  (27))  become  applicable  for  pa  greater  than  5). 

Hand  computations  for  a  few  specific  instances  demonstrate  this  phe¬ 
nomenon. 

Example  1.  For  a  ceramic  material  1/4  cm  thick  in  turbulent  flow,  under 
a  beam  of  10  cm  diameter,  taking  h  =  10~3  watt/cm*°C  and  k  =  -.01  Joule/ 
sec  cm°C,  we  find  pa  “  6.  In  this  case,  r*/a  “  1.1  and  stresses  are  found 
as  follows  from  Figure  1  and  Equation  (31). 


ils. 

EE 

1.12 

6.8 

-0 

1.16 

7 

.093 

1.2 

7.2 

.145 

1.4 

8.4 

.257 

1.5 

9 

.256 

1.6 

9.6 

.231 

2.0 

12 

.169 

Example  2.  A  disk  of  alumina  (AlgOj)  1/2  cm  thick  in  the  same  turbulent  flow 
(h  =  10  watt/cm*  C)  with  k  -  .2  Joule/sec  cm  °C  and  a  6  cm  diameter  beam 
leads  to  pa  =  .4.  In  this  case,  tensile  stresses  develop  at  pr*  =  1.3.  A 
few  stresses  are 
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lL± 

RE 

-3 

1.3 

-0 

3.5 

1.4 

.0161 

5 

2.0 

.0507 

6.25 

2.5 

.0547 

7.5 

3.0 

.0519 

9.0 

3.6 

.045 

Example  3.  A  disk  of  aluminum  (k  =  2  Joule/sec  cm  °C) ,  otherwise  identical  to 
the  preceding  example  leads  to  pa  ~  .1.  In  this  case,  tensile  stresses  de¬ 
velop  for  pr*  =  1.2  or  r/a  =  12.  Some  stresses  are 


r/a 

£1 

°ee/EoT 

12 

1.2 

-0 

14 

1.4 

.0135 

20 

2.0 

.0267 

22 

2.2 

.0276 

25 

2.5 

.0275 

30 

3.0 

.0257 

50 

5.0 

.0144 

For  very  good  conductors  k  *  00  and  pa  0.  In  this  case  pr*  -*■  1.2  and 
r*  *  00  as  k  ♦  00.  Since  Lim^^g  F^(pa)  1,  we  may  develop  a  simple  bound 
on  the  maximum  tensile  stress  possible  in  a  good  conductor.  Setting  KQ(p«)  = 
-ln(pa)  and 


FL(pa)  =  1 

in  Equation  (31),  and  differentiating  with  respect  to  r,  we  find  that  the 
maximum  tensile  stress  always  occurs  at  pr  =  2.32.  Hence,  a  bound  on  the 
tensile  stress,  valid  for  pa  <  .3,  is 


=  _  0.0687 
EoTj^  "  ln(pa) 

max 


(32) 
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Thus  the  tensile  stress  always  develops  but  at  Increasingly  large  (r/a  = 


2.32/pa)  distances  fro*  the  heated  area. 

It  Is  easy  to  show  that  for  large  pa  (the  perfect  insulator)  the  maximum 
tensile  stress  occurs  just  outside  the  heated  area  and  is  of  magnitude 


These  computed  values,  together  with  the  stress  distribution  for  a  per¬ 


fect  insulator  (k  =  0)  are  sketched  in  Figure  2.  The  decrease  in  the  maximum 
tensile  stress  for  decreasing  pa  (increasing  conductivity)  is  evident.  The 
bound  (Equation  (32))  is  given  as  the  dashed  line  and  represents  the  envelope 
enclosing  all  peaks  obtained  for  pa  <  0.3. 


Figure  2.  Stress  Distribution  (c  )  for  Various  Values  of  p 


B.  Quasi-Static  Thermal  Stresses 


1.  Local  Heating  of  Thin  Sheets 

As  a  first  example  of  a  quasi-static  thermal  stress  distribution,  let  us 
consider  the  process  whereby  the  long-time  stress  distribution  given  in  the 
preceding  section  is  developed.  We  first  consider  the  special  case  of  the 
poor  conductor,  for  which  the  temperature  outside  of  the  heated  area  remains 
at  the  ambient  temperature,  but  assume  a  uniform  temperature  T^(t)  for  r  <  a. 
These  assumptions  are  appropriate  for  heating  times  short  compared  to  radial 
diffusion  times, 

to  <<C  fcd  =  a*/k  (1) 

but  long  compared  to  the  thermal  diffusion  time  through  a  thin  sheet 

tQ  «  l*/k  (2) 

In  this  case 

T»  -  -  js*r  ‘  <3> 

p 

for  a  beam  of  uniform  absorbed  intensity  1^. 

In  consequence  of  this  temperature  distribution,  thermal  expansion  occurs 

ever  o  <  r  <  a  which  is  partially  restrained  by  the  surrounding  material  which 

in  still  (since  we  are  ignoring  thermal  conduction)  at  T  =  T  .  Thus,  the 

o 

stress  Held  which  results  is  that  resulting  from  forcing  a  disk  of  radius 

a’  =  a(  1  +  oT^  (4) 

into  a  hole  of  radius  a  in  a  sheet  of  the  same  thickness.  The  displacements 
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in  the  inner  and  outer  regions  are  the  well  known  radial  distributions  of  the 
Lame  problem* 


U  =  Ar  +  B/r 

(5a) 

U  =  Cr  +  D/r 

(5b) 

0 

matching  normal  stress  on  r  =  a  and  requiring  continuity 

a '  +  U .  =  a  +  U 

(6) 

i  o 

leads  to  C  =  B  =  0  and 

D  =  a*a  (-^y-*)  T, 

(7) 

A  =  -a  ( 1  2  V)Tl 

(8) 

From  which  we  find 

EaT^ 

(9) 

arri  =  '  “ 

EoTi  a* 

°rr  =  "  —  (7r) 

(10) 

0 

In  addition  to  these  radial  stresses,  circumferential  stresses 

also  develop 

and  are  of  magnitude 

EaT^ 

a66i  =  2~ 

(ID 

EaT)  .* 

°60  "  2  <S*> 

0 

(12) 
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for  the  assumed  state  of  plane  stress.  It  Is  these  tensile  stresses  outside 
of  the  heated  region  which  we  may  expect  to  have  an  effect  on  crack  propaga¬ 
tion.  The  maximum  value  is 


EaTj^ 

aee  =  2 

max 


(13) 


In  a  ductile  material,  some  yielding  is  to  be  expected  when  the  maximum  shear 
stress  exceeds  a  critical  value.  In  the  heated  region, 


EaT^ 

Tmax^  -  4 

everywhere,  and  outside  the  heated  region 

EaTL 

Tmax  ~  2 

o 


(14) 


(15) 


at  the  edge  of  the  heated  zone. 

We  may  note  a  fracture  criterion  for  a  brittle  material  which  is  a  poor 
thermal  conductor.  Under  these  conditions,  a  uniform  deposition  of  energy  Q 
over  a  circle  of  radius  a  will  lead  to  fracture  if 


pC  1 

q  _  2o  a*  ( — 2--) 

xrlt  crita  v  Ea  ' 


(16) 


and  °crit  is  the  tensile  stress  producing  fracture.- 

We  may  tak?  thermal  softening  into  account.  If  the  modulus  of  the  heated 
material  is 


E(T1)  =  E1 


(17) 


22 


23 


Thus,  while  therm* 1  softening  relieves  some  of  the  stresses,  it  does  not 
eliminate  them  coaqiletely. 

As  a  demonstration  of  the  manner  in  which  thermal  softening  within  the 
heated  region  curtails  the  development  of  tensile  stresses  outside,  let  us 
consider  a  hypothetical  material,  for  which 

E(T)  =  Eocos(|  J-)  (25) 

m 

Heating  Is  presumed  to  occur  in  accord  with  Equation  (3),  and  tensile  stresses 
as  given  by  Equation  (22)  to  develop. 

We  find 


90 


max 


EoT 


=  —  [(1  +  v)  +  (1  -  v)  sec  (~  —)] 

m  m 


-1 


(26) 


where 


t  =  PC  IT  /I  (27) 

n  p  b  a 

is  the  time  required  to  produce  meeting.  This  result  is  depicted  graphically 
as  Figure  3, 

Let  us  now  consider  the  occurrence  of  yielding  in  the  heated  region  of 
an  elastic,  perfectly  plastic  material  of  negligible  thermal  conductivity, 
recognising  that  these  two  assumptions  are  in  large  part  mutually  exclusive, 
for  those  materials  which  yield  plastically,  e.g.  metals,  tend  to  be  good 
thermal  conductors.  As  is  well  known,  ductile  metals  have  yield  points  which 
are  strongly  temperature  dependent.  Thus,  if  the  process  described  in  the 
preceding  progt esses  to  some  temperature  such  that 
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WITHOUT  SOFTENING 
WITH  SOFTENING 


Figure  3.  Thermal  Stress  (Tensile)  in  Hypothetical  Material 


a  a 


i 


(28) 


EaT. 


=  o  (T  ) 

y  1 


further  heating  will  not  increase  the  stresses  for  either  r  <  a  or  r  >  a. 

Rather,  further  heating  will  reduce  the  yield  points  and  hence  the  stresses. 

If  heating  is  stopped  at  some  temperature  T,  (T,  <  T  <  T  ),  the  stress  dis- 

xi  i  in 

tribution  at  that  time  will  be 


-VV 


(29) 


if  thermal  softening  and  heat  conduction  for  r  >  a  is  neglected.  We  now  sup¬ 
pose  that  the  entire  plate  returns  to  the  ambient  temperature,  and  a  process 
of  elastic  recovery  occurs.  The  stresses  which  develop  during  this  process 
are  readily  seen  to  be 


-  Ea(Tf )/2 


a* /r* ) 


(30) 

(31) 


superposing  these  onto  the  stresses  at  the  end  of  the  heating  period  gives  the 
final  state  of  stress  inside  r  =  a  as 


Ea  -f  -  o  <T,) 
2  y  f 


and  for  r  =  a  in  the  exterior 


(32) 


f  -  ,y<T{> 


(13) 


Thus,  it  is  seen  that  large  residual  tensile  stresses  can  develop.  If 
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the  elastic  perfectly  plastic  material  is  heated  to  the  melting  point,  and  the 
melt  remains  in  place,  these  final  stresses  can  be  as  large  as 


a 

max 


=  Eo  T  /2 

ID 


(34) 


which  may  be  of  the  order  of  the  yield  stress  at  room  temperature.  Moreover, 
these  large  tensile  residual  stresses  may  produce  fracture.  Thus,  a  single 
large  pulse,  not  quite  capable  of  producing  melt  through,  may  give  rise  to 
mechanical  failure.  The  time  history  of  the  process  is  depicted  graphically 
in  Figure  4.  The  process  occurs  as  follows.  During  0  <  t  <  t^,  heating 
occurs,  (OA)  and  elastic  stresses  equal  to  o^(T^)  develop.  From  t^  <  t  <  t^, 
heating  continues,  but  the  heated  disk  deforms  plastically  (AB)  and  the 
stresses  fall  to  o^(T?).  For  t  >  t^,  cooling  takes  place,  and  the  material 
loads  elastically  (BC)  to  the  final  stress,  given  by  Equation  (33).  In  a 
work  hardening  material,  the  stress  at  t^  will  be  greater  and  the  final  value, 
of,  reduced. 

This  damage  mechanism  appears  to  be  of  such  significance  as  to  merit 
further  study.  In  particular,  the  degree  to  which  radial  conduction  away 
from  the  heated  region  and  strain  hardening  alleviates  the  stresses  should 
be  examined. 

Proceeding  in  the  same  manner  as  in  Section  A,  it  can  be  shown  that  con¬ 
duction  for  r  >  a  during  heating  does  not  affect  the  state  of  residual  stress 
after  cesssation  of  heating  and  the  establishment  of  thermal  equilibrium. 

The  influence  of  strain  hardening  may  be  developed  as  follows.  We  assume 
the  stress  strain  curve  for  the  material  is 

o  =  Ee  for  o  <  a  (35a) 

y 

a=o  +  8E(e  -  c  /E)  for  a  >  a  (35b) 

y  y  y 
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time 


where  is  taken  to  be  a  function  of  temperature  but  E  and  0  are  not.  This 
stress-strain  curve  is  depicted  graphically  as  Figure  5. 

Let  us  again  suppose  the  region  r  <  a  is  again  raised  to  a  temperature 
T^  and  that  yield  initiates  uniformly  on  r  <  a  according  to  a  maximum  shear 
criterion.  He  then  assume  that  heating  continues  to  a  temperature  T2  >  T ^ 
with  the  final  strains  being 


E  =  E 
rr  ee 


Hence 


°oei  =  °rri  =  °y(T2)(1  “  6)  +  BE  a 


These  radial  stresses  and  displacements  must  match  (at  r  =  a)  those  of  the 
exterior  region,  =  -C^E/[( l+v)r* ]  and  Ur  =  C^/r  at  r  =  a.  Hence 


+o  (T2)(l  -  8)a* 


4  E(-B  +  •—> 
1+v 


and  the  stresses  in  r  <  a  become 


rt  =  °«8j  =  [l~  (1  + 


-oy(T2)(l  -  8) 
[1  -  (1  +  v)B 1 


while  for  r  >  a 


-o  (T2)(l  -  8)  >t 


rro  [T  -  8(1  +  v)]  r1 


(40a) 


oy(T2)(l  -  8)  a, 
“  [l  -  8(1  +  v)]  r* 


(40b) 
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Figure  5.  Stress-Strain  Relationship  for  Hypothetical  Strain-Hardening 
Material . 


Cooling  to  the  ambient  temperature  again  leads  to  the  superposition  of  Lame 
stresses,  and  the  final  residual  stress  state  is: 


Ea(T2)  ov<T2)(1  -  8) 
°rri  =  °ee1  =  2  [1  -  (fTT)H 


(41a) 


EaT2  o  (T2)(l  -  B)  #t 
°rr  =  “°ee  =  {“T"  “  [¥-  (i  +“)b]>  7* 


(41b) 


We  conclude  that  strain  hardening  will  reduce  the  magnitude  of  the  reside 
ual  stresses  appreciably  only  if  the  final  temperature,  T2,  is  only  slightly 
above  the  value  at  which  yielding  occurred.  For  example,  if  the  relationship 
between  yield  stress  is  taken  to  be 


o  =  a  (1 - ^ - ) 

y  yo  1200°C 

(42) 

which,  at  least  qualitatively,  describes  data  available* for 

several  materials. 

V  =  i 

(43) 

(44) 

the  residual  stress  without  strain  hardening  is 

E°T2  „  T2  , 
oB=  „  -a  ( 1  -  ) 

R  2  yo  1200  C 

(45) 

*Holmes,  fi*  S.,  and  Desmond,  T.  P.,  Thermal  Mechanical  Damage  Study, 
Monthly  Status  Report  #8,  Contract  #  F29601-78-C-0041 ,  SRI  International, 
May,  1979. 
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and 


If  strain  hardening  la  Included. 


2.  Uniform  Heating  of  Platea 

In  the  preceedlng  section,  thermal  stresses  resulting 
from  radial  variations  in  temperature  were  considered,  while  a 
uniform  temperature  was  assumed  through  the  thickness.  We  turn 
now  to  the  other  limiting  case  -  that  of  variations  in  temperatures 
through  the  thickness  of  a  plate,  with  no  temperature  variations 
in  the  transverse  direction,  as  would  be  encountered  if  energy  is 
absorbed  at  a  uniform  rate,  Ia,  over  a  sufficiently  large  area, 
a.  The  Free  Plate. 

We  first  assume  the  plate  to  be  free,  that  is,  that 
there  are  no  boundary  constraints.  For  convenience,  we  take  the 
plate  to  be  rectangular  of  thickness  h,  bounded  by  planes  Z**0  and 
Z-h. 


Ignoring  mechanical 
plane  stress  will  result,  i.e,  0 
symmetry,  we  see  also  that  ©xy  ■ 
stress  strain  relationships  then 


inertia,  we  see  that  a  state  of 

2  ■  0.  Th'.ough  considerations  of 

0  and  that  0  -  0  -  o(Z).  The 

xx  yy 

give 


Eexx  *  0-'v)  a<Z>  +  T  -  EF (Z)  (1) 

Eeyy  -  (1-v)  0 (Z )  +  Ea  T  «  EF(Z)  (2) 

and  inplane  displacements  may  be  obtained  through  integration  of 
the  strains.  The  result  is  of  the  form: 

Ux  -  x  F (Z)  +  h(Z,  y)  (3) 

Uy  -  y  F(Z)  +  g ( Z ,  x)  (A) 

The  thickness  strain,  and  displacement  are  (respectively) 

eZZ  “  "  ^E  °(Z)  +  a  T  *  G(Z)  (5) 

U2  -  /G(Z)dZ  +  k(x ,  y)  (6) 

Shear  stresses  are  then  seen  to  be 

-  P  U  P'(Z)  +&♦!£>  (7) 

V  -  p  <y  f’U)  +  |f  +  I?1  (8) 
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Since  o  *  a  ■  o(Z)  and  o  •  0,  Che  equilibrium  equation! 
xx  yy  xy  ’  n  ’ 

require  that 


3o 


—  -  0  or  XF"(Z)  +  -  0 

az2 


3  Z 

3o  z  32g 

— *=■  -  0  or  yF"(Z)  +  -  -  0 

3Z  3x2 

Since  h  is  independent  of  x,  and  g  of  y,  we  see  that 
F’  »  const  ■  Cj.  Hence 

F(Z)  -  Ci  Z  +  C2 

and  we  arrive  at  a  stress  dlsbribution 


(9) 


(10) 


(ID 


o(Z)«  a  -  cr  -  {Ci  Z  +  C2  -  aT}  (12) 

xx  yy  1-v  1  *■ 

where  temperature,  T,  is  measured  with  respect  to  an  ambient, 
unstressed,  state. 

If  we  require  the  plate  edges  to  be  free,  the  moment  per 
unit  length  muse  vanish,  i.e., 

MT  -  0  -  /h  o(Z  -  h/2)  dZ  -  0  (13) 

o 

and  the  in-plane  force  per  unit  length  must  also  vanish,  i.e., 

Ft  -  /h  adZ  -  0  (14) 

1  o 

These  give 


C2 


—  /h  (2h  -  3Z)  TdZ 
h2  ° 


(15a) 


Cl 


—  /h  (2Z  -  h)  TdZ 
h3  ° 


(15b) 


From  which  we  may  deduce  the  final  expression  for  the  thermal  stress 
in  a  free  plate,  uniformly  heated  over  one  entire  face,  tc  be 

o  -  o  -  f2-|  —  /h  (2z-h)  TdZ  +  —  /“  (2h  -  3Z)  TdZ  -  t"|  (16) 

xx  yy  0  h*  °  J 


It  is  readily  shown  that  a  uniform  or  linear  temperature  produces 
no  stress,  even  if  the  coefficients  are  time  dependent,  as  in 

T  -  A(t)Z  +  B (t )  (17) 
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For  a  plate  absorbing  energy  at  a  rate  Ia  per  unit  area, 
per  unit  time,  on  the  plane  Z  -  h,  with  no  losses  on  the  f;  :e 
Z  -  0,  the  temperature  may  be  shown  to  be 


It  I  h 

T  (Z  ,  t )  -  — -  +  -5-  (Z2/h2-l/3) 

pCph  2k 


£ 

n-1 


21  h 
a 

kn2ir2 


,  .  v  n  -icn2Tr2t/h2 
V-l )  e 


For  large  values  of  time,  the  series  makes  negligible  contributions, 
and  the  temperature  may  be  satisfactorily  approximated1  by  the  first 
two  terms.  Of  these,  the  first  is  recognized  as  leading  to  no  stress. 
Thus,  for  long  times,  the  stresses  are  generated  by  the  quadzatic 
term  in  the  temperature,  and  are 

Eu  I  h 

a  «  a  ■  - — — 

XX  yy  2  (1-v) k 

The  maximum  value  occurs  at  the  front  and  rear  faces,  and  is  a 
compressive  stress  of  magnitude 

Ea  I  h 

a  - - - -  (20) 

12 (1-v) k 


Z 

h 


Z2 

h2 


(19) 


The  maximum  tensile  stress  occurs  at  the  mid-plane,  and  is  of 
magnitude 


Ea  I  h 

a  -  - - -  (21) 

24 (1-v) k 

These,  we  note,  are  static  stresses  even  though  the  generating 
temperature  (Equation  18)  increases  linearly,  with  time,  at  long 
times . 


1 

Torvik,  F.J.,  Thermal  Response  Calculations  and  Their  Role  in  the 
Design  of  Experiments,  AFIT-TR.-7 3-6  .  Also  in  Proceedings  of 
1973  DOD  Laser  Effects  and  Hardening  Conference,  Monterey  CA, 

Oct  1973. 
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rinyfnl  - 


>- 


At  short  times,  stresses  generated  by  the  series  terms  must  also  be 
Included 


0trans 


E®  bfz  „(-l)n  -ten *  it*  t  / h  *  nir  Z  /2Z  N24  „  1  -icn2n2t/h2 

— r — p— £-*  f  e  cos  —  -i  —  -  IF-  I  —  e 

k(!-v)  p  q2  h  l^h  ,Vodd  n4 


This  result,  with  some  notatlonel  changes.  Is  to  be  found  In  Boley1 
and  has  been  shown  to  be  negligible  for  slowly  varying  heat  Inputs. 
In  heat  additions  which  are  as  step  functions  In  time,  as  are  of 
Interest  here,  these  terms  are  of  significance  for  short  times. 

Some  numerical  results,  obtained  by  Garrison2,  are  presented 
in  Figures  6  and  7.  In  figure  6  the  distribution  of  stress  through 
the  thickness  is  shown  for  several  instants  of  time.  The  general 
increase  with  time  of  the  peak  tensile  stress  and  the  movement  of 
the  place  of  occurence  to  the  center  of  the  plate  can  be  seen.  At 
Kt/h2  >  1,  the  stresses  are  Indistinguishable  from  the  static  value. 
The  compressive  stress  at  the  heated  surface  first  Increased  and 
then  decreased.  These  temporal  changes  are  more  clearly  seen  in 
Figure  7.  The  peak  compressive  stress  is  seen  to  reach  a  value 
in  a  finite  plate  some  30%  greater  than  the  asymptotic  value  for 
long  times,  as  given  by  Equation  19. 

b.  The  Influence  of  Constraints. 

This  simple  example  of  one  dimensional  temperature  variation 
in  a  plate  presents  an  opportunity  to  investigate  the  influence  of 
mechanical  constants  at  the  boundary  on  the  stress  distribution. 
Computing  the  solutions  for  the  displacement  (using  the  symmetry  of 
x  and  y  axes),  we  find 


22 


Ux 

uy 


Ci  x  Z  +  C  2X  4-  D 
CiyZ  +  Cjy  +  D 


- 

;o  1- 


(l  +  v) 


aT  -  —■ 

\>  1- v 


(2  3a) 
(23b) 


C 

(Ci Z+C 2)  )  dz  -  j*-  (x2+  y2)  +  E  (23c) 


1Boley  and  Wiener,  p  285. 

2Garrison,  Jan  N.,  Thermal  Stresses  as  a  Laser  Heating  Damage 
Mechanism,  AFIT  Thesis,  GAE/MC/75-4,  November  1976. 
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COMPRESSIVE  STRESS 


ik  Thermal  Stress  Versus  Dimensionless  Heating  Time  for 
!— Dimensional  Axial  Heat  Flux  (Dimensionless). 


As  our  example,  we  consider  a  square  plate  with  side  length 
a,  bounded  by  x  ■  0,  y  *  0,  x  ■  a,  y  ■  a.  Pour  possible  boundary 
conditions  are  possible: 

I  The  edges  free. 

II  The  edges  constrained  to  prevent  rotation,  but  in-plane 
expansions  freely  allowed. 

Ill  The  edges  constrained  against  ln-plane  displacement  at 
the  mid-plane,  but  rotations  allowed. 

IV  The  edges  constrained  against  both  displacement  and 
rotations . 

The  stresses  for  case  I  are  those  of  Equation  16. 

If  no  rotations  are  allowed,  inspection  of  the  displacements 
reveal  that  Ci  -  0.  Ensuring  the  absence  of  in-plane  forces 
requires 

/h  odZ  -  0  (24) 

o 

or 

C2  -  ~  /  TdZ  (24) 

hence  the  stresses  are,  for  case  II, 


a  -  a  -  /h  TdZ  -  T) 

xx  yy  1  —  v  h  o 


(25) 


If  now  we  require  that  <n-plane  displacements  vanish,  for 
Z  ■  h/2,  on  the  boundary,  we  see.  iron  equation  23  that 

C2  -  -  \  h/2  (26) 

Vanishing  of  edge  moment  requires  th;  t 

/h  (Z  -  h/l>  a  (4)  dZ  -  0  (13) 

o 

or 

Cj  -  -  /h  TdZ  -  fh  ZTdZ)  (27a) 

3  '  ^  O  O 

C 2  m  —  (%  /h  TdZ  -  fh  ZTdZ)  (27b) 

2  \  2  o  o 

h 

Finally,  for  case  III. 
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0__  "  “  f^r  TdZ  ‘  /JjZTdzV  —  /h  TdZ  -  ;h  ZTdz)  -  T }  (28 ) 

x*  yy  l-v  o  o  /  ^3  \2  o  o  / 


The  fourth  case  arlsaa  for  0X  »  By  ■  0  for  all  x,  y,  t,  or 
Cj  -  0"  C2  ■  0.  Than 


0 


xx 


0 

yy 


Ea 

l-v 


<-T) 


In  ordar  to  aaseaa  the  Influence  of  edge  constraints  on  the  stresses, 
let  us  put  the  temperature  distribution  for  long  tines  (Equations  18) 
into  the  forn 

T(Z , t )  -  ~y~  Q  +  7  (z2/h2  “  i/s)3  (29) 

where 


Kt 


is  the  dimensionless  tine  previously  introduced 
Then  the  stresses  are  found  to  be,  in  each  case: 

I.  For  no  edge  constraints: 


(31) 


II.  For  edges  constrained  to  prohibit  rotation  only  (expansion 
permitted) 


o 


XX 


(32) 


III.  For  edges  constrained  to  prohibit  expansion  only  (rotation 
permitted) 

0  ■  o  -  T^hv-r-)  (-  t  -  —  +  1/2  Z/h  -  1/12}  (33) 

xx  yy  i-v\  k  /  2h2 

IV.  For  edges  constrained  to  prohibit  both  rotation  and  displace- 
mei’ 


xx 


yy 


„  Ea  /XahV  .1  Z2 
"  ){-  T  +  6  -  —2 


(34) 
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In  all  cases  the  maximum  stress  Is  compressive,  and  occurs  at 
the  heated  face  and  Is  Independent  of  plate  site.  We  note  that 
it  is  the  in-plane  constraints  rather  than  the  bending  constraint 
which  generates  the  large  stresses  at  long  times.  At  shorter  times 
however,  t  'v  .25  or  less,  we  note  that  greater  stresses  are 
generated  through  the  bending  constraint.  Although  these  simple 
expressions  are  not  valid  for  values  of  t  much  smaller  than  this 
value,  the  same  generalization  is  valid  for  very  small  t. 

At  very  short  times,  t  <  <  1,  we  may  regard  the  plate  as 
being  of  infinite  thickness.  In  this  case, 

T  (Z  ,  t )  -  —5—  /Tt  ierfc  (- - ^  (35) 

k  \2/<t/ 

where  Z  is  now  measured  into  the  slab>  with  Z  ■  0  being  the  heated 
face.  Recognizing  that  the  unheated  (and  therefore  unexpanded) 
mass  of  thick  plate  serves  as  a  constraint  on  the  heated  surface 
layer,  the  stresses  are  (from  case  IV) 


a 

xx 


a 

yy 


EaT 

1-v 


(36) 


The  maximum  stress  is  again  compressive,  at  the  surface,  and  given 
by : 


a  (0,t)  -  0  (0,t) 

xx  *  yy 


Ea  2Ia 
1-v  k 


(37) 


The  dashed  line  of  Figure  6  shows  the  stress  at  the  surface 
computed  from  Equation  37,  These  values,  which  are  analogous  to 
complete  constraint,  are  slightly  greater,  at  small  times  than  the 
numerical  results  obtained  for  the  free  plate.  The  stresses  in 
the  thick  slab  decay  rapidly  with  depth  according  to 


0  «o 
xx  yy 


Ea 

1-v 


21 


/IcT  ierfc 


(38) 


with  values  at  times  comparable  to  «tct/h2  ■  .0057.  Computed 
results  have  been  added  (in  dashed  lines)  to  Figure  7  for  purposes 
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of  comparison  with  ths  rssults  obtained  for  tha  plate  of  finite 
thickness.  The  rapid  decay  from  the  surface  is  evident.  From 
a  table  of  values  of  the  ierfc  function,  it  nay  be  deduced  that 
the  stresses  fall  to  10X  of  the  surface  value  at  a  depth  Z 
'v  2/ITt,  and  to  IX  of  the  surface  value  at  a  depth  of  Z  ^  3/<"t. 


C .  Thermal  Shook 


As  a  first  example  of  a  problem  wherein  the  effects  of  thermal  inertia, 
as  well  as  mechanical  inertia,  are  significant,  we  consider  a  thin  elastic 
plate,  circular  in  planform,  and  supported  at  the  outer  periphery  (clamped, 
or  simply  supported).  A  uniform  energy  flux  I  is  absorbed  over  the  entire 
surface,  commencing  at  t  =  0.  Thus  we  assume  a  through-the-thickness  temp¬ 
erature  variation  only,  which  generates  a  thermal  moment.  The  response  of 
the  plate  to  the  uniform,  but  time  varying  moment  is  desired. 

We  let  w(r,t)  be  the  deflection  of  the  mid  plane  of  a  plate  bounded 
by  -h/2  <  z  <  h.  Then,  using  elementary  plate  theory 

ur  =  "  z  fr  5  ue  =  0  (D 

3*w  z  aw 

crr  Z  aY*  ’  cse  ”  r  ar 

Hence,  for  a  linear  elastic  material  at  temperature  T  above  the  ambient 
(unstressed)  temperature 


-  ~E  '  z 

3*W  vz  3w  1  E  _ 

/  1  Q  \ 

rr 

i-vT 

ar*  +  r  8r  J  1-v 

V  *5®  ) 

E  r 

£  12  +  i!"  1  _  JL  aT 

i 

ree 

1 

I 

»-» 

1 

< 

r  ar  ar*  J  1-v 

We  define  total  moments  (per  unit  length)  and  a  thermal  moment,  according  to 


fh/2 

«r  =  f  J  zorrdz 
-h/2 

fh/2 

"«  ’  %edz 


zEaTdz 


-h/2 


(4a) 

(4b) 

(4c) 


43 


This  produces 


N 


r 


-  D{ 


>*W  v  9w 
•r*  +  r  9r‘ 


"T 

1=7 


M 


e 


D 


+v 

r  9r 


(5s) 

(5b) 


where 


D 


Eh* 

l2(  W) 


(8) 


An  obvious  division  of  the  total  moments  into  mechanical  end  thermal  compo¬ 
nents  can  be  made.  Since  no  elongation  of  fibers  on  the  neutral  plane  of 
symmetry  is  allowed  in  elementary  theory,  in-plane  stresses  will  also  develop. 

From  the  sign  conventions  and  free  body  diagrams  of  Figure  8  we  deduce 
that  equilibrium  of  moments  about  the  0  axis  requires  that: 

(N  +  dr)(r  +  dr)de  -  M  rde  -  (M  +  ~S  rde)drde  -  Qrdrde  =  0  (7) 

r  #r  r  0  8“ 

or 

r  +  (Mr  -  Mft)  -  Qr  »  0  (8) 

Summing  vertical  forces  gives  (from  Fig  Id) 

fr 

I  ph2ww(r,t)rdr=  2wrQ  (9) 

o 

Thus  one  form  of  the  governing  equation  for  the  plate  is: 

3N  M  -  M  .  ,r 

— 2_)  -  -  j  Phrw(r,t)dr  =  0  (10) 

o 

We  note,  from  5,  that 

»r  - ",  *  -  D(1-',)r  a?  <?  <111 


We  note  differentiation  produces 

,r 


~  |  phrw(r,t)dr  =  phrW(r,t) 


(12) 
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Figure  8.  Plate  Elements. 
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and  that 
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has  the  solution 


a*w  v  aw 

l?  +  ?Fs0  on  r  =  a 

(20a) 

WH  =  0  on  r  =  a 

(20b) 

wH(r,o)  =  f (r) 

(21a) 

wM(r,o)  =  g(r) 

(21b) 

wu  =  y  {A  sin  to  t  +  B  cos  u  t}R  (r) 
n  n  n  n  n  n 

n=l 


where  the  eigenfunctions  R^fr)  are 


J  (p  a) 
o  n 


with 


Rn(r)  =  Jo(Pnr)  “  nr^T  Io(Pnr) 
o  n 

fV  _  Jl(pnal 

1-v  J  (pj*a)  +  I  (p*  a) 
o  n  o  n 


The  coefficients  A^,  are  given  by 


A  to 
n  n 

-  i-  fr 

g(r) 

B 

n 

“  N  P 

n  o 

f(r) 

R  (r)dr 
n 


(22) 


(23) 


(24) 


(25) 


where  Nr  -  j  r{Rn(r)}2dr  (26) 

o 


Returning  to  Equation  16,  we  expand  the  left  hand  side  in  terms  of  the 
eigenfunctions : 


Let 

tb  ’2»t  ■  2, 

(27) 

n-1 

and  let 

w(r,t)  -  w*(r,t)  +  p(r)q(t) 

(28) 

where  the  functions  p(r)  and  q(t)  are  to  be  chosen  so  as  to  leave  homogeneous 
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, ,  .  „„  „*  at  r  -  a.  Such  a  choice  is 

boundary  conditions  on 

p(r)  -  (a  -  r) 

(29) 

a  M  (a,t) 

q(t) - f - 

\>D  1-v 

(30) 

Let 

ph(a-r)q  -  I  G^OR^V5 
n«l 

(31) 

and 

w*(r,t)  -  I  Wn(t)Rn(pnr) 

(32) 

Then  lubstitution  yield* 


phWn(t)  +  -  -  (F„  +  °n> 

(33) 

where 

p  <t)  .-i-i-  [*t^HT<I-t)Rn<V)dI 

nvw  w  1-v  J  1 

n  o 

(34) 

and 

G  <t)  .  _1  [‘r(»-ri*„(V>ar 

GnW  Nn  vD  l-°  £ 

(35) 

The  new 

initial  conditions  become: 

w(r,o)  -  0  -  w*(r,o)  f  (a-r)q(o) 

;(r,0)  -  0  -  w*(r,o)  +  (a-r)q(o) 

(36a) 

(36b) 

Thus 

a  Mj^a.o) 

I  wn<0)Rn(pnr)  “  "  (a"r)  >®  1_v 

(37a) 

a  M^(a,o)_ 

l  V0)  '  Rn(pnr)  "  ’  (a”r)  "  1"v 

(37b) 

or 

„  (0)  .  -  i-  i-  f  rC-r>*n<»»r)4t 

Wn'0J  vD  1-v  N  J 

n  0 

(38a) 

0 

(38b) 
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Given  a  time  history  of  temperature,  then,  may  be  computed  from  Equation 
4c,  and  the  time  dependent  amplitudes  determined  from  equation  33,  34  and  35, 
with  initial  conditions  given  by  38. 

The  solution  for  the  analogous  problem  of  the  beam  and  the  rectangular 
plate  have  been  given  by  Boley*  and  the  observation  made  that  the  results 
for  the  thermal  vibrations  of  beams  differ  but  little  from  those  of  plates. 
That  being  the  case,  it  would  appear  that  the  beam  results  should  suffice 
for  the  circular  plate  as  well,  if  the  effective  beam  length  is  taken  to  be 
diameter.  The  results  are  expressible  in  terms  of  a  ratio  of  dynamic  to 
static  deflection, 


1 

w  _  1  +  2e  for  B.  >  -  .69 

-  1  - 

Wmax  static  2  for  B^  <  -  .69 


(39) 


where 


h 

2aV*T 


<*-> 

hp 


1/4 


(40) 


and  a  is  the  plate  radius.  Physically,  Bj2  corresponds  to  the  ratio  of  thermal 
response  time  to  mechanical  response  time.  The  static  deflection  of  a  simply 
supported  plate  under  thermal  load  is  obtained  from  Equation  16: 


DV^w  +  —  *  harmonic  function  of  r  only  (41) 

1-v 

For  a  solid  plate,  with  temperature  Independent  of  r,  and  having  a  3imply 
supported  edge. 


w 


MT(a2  -  r2)6 
Eh3 


The  stresses  are 

12M^,  E 

o  ■  ~7Z — r— r  i  -  - —  oT 
rr  (l-v)h3  1-v 


1 

Boley  and  Weiner,  pp  406-409. 


(42a) 


.■IT-C.-tWMl 
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(42b) 


_  12MT2  EaT 

°ee  “  (l-v)h*  -  1-v 

These,  it  should  be  noted,  are  Identical  to  the  stresses  obtained  in 
Case  III  of  section  B2b,  aside  from  changes  in  notation. 

To  account  for  the  dynamic  effect  (approximately),  we  may  compute  the 
dimensionless  ratio  B^,  from  equation  40,  and  the  stresses  from  42,  and  adjust 
by  multiplying  by  the  ratio  R  computed  from  equation  39. 

In  the  case  of  clamped  edges,  we  note  that  the  boundary  conditions  are 
naturally  homogeneous,  even  in  the  presence  of  the  thermal  loading,  i.e., 

w(a,t)  =  g  (a,t)  -  0  . 

Hence,  for  homogeneous  initial  conditions,  we  need  only  solve  equation  16. 

We  note  immediately  that  in  the  special  case  where  the  temperature  is  inde- 
oendent  of  r,  that 

=  0 

Hence,  w(r,t)  =  0  and  the  stresses  for  the  dynamic  problem  (thermal  shock) 
are  the  same  as  for  the  quasistatic  problem,  as  no  displacements  develop. 

This  is  correct  only  to  the  limits  of  thin  plate  theory,  which  does  not  pro¬ 
vide  for  "through  the  thickness  displacements"  in  addition  to  those  of  the 


neutral  axis. 


III.  STRESS  GENERATION  BY  AIR  ABSORPTION 


A.  The  Generation  of  Front  Surface  Pressures 

We  consider  a  laser  beam  absorbed  over  a  surface  of  a  progressive  front , 
leaving  behind  a  pressure  P,  density  p,  particle  velocity  U,  measured  in  an 
inertial  frame.  The  front  is  taken  to  propagate  at  speed  D.  By  conservation 
of  mass 


PqD  =  pf(D-U> 


By  conservation  of  momentum 


?r  -P  = 
f  o 


p  DU 
o 


(1) 


(2) 


By  conservation  of  energy 


AE 


P 

o 


(3) 


Aside  from  AE,  the  energy  added  per  unit  mass,  these  are  the  familiar  Rankine 
Hugoniot  Jump  conditions  across  a  shock  front. 

Here,  the  added  energy  per  unit  mass  is  computed  by  observing  that,  for 
any  area,  A,  of  the  front,  the  laser  flux  f  arriving  over  time  At  is  assumed 
to  be  completely  absorbed  by  the  volume  of  material  processed  during  the  same 
Interval,  i.e. 


AE  •  fA  At/Am 


The  mass  absorbing  AE  is 


Am  =  Ap  DAt, 
o 


(4) 


(5) 
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Hence 


AE  =  f/p  D 
o 


(6) 


For  a  perfect  gas. 


E  =  CVT 


P  =  pRT 


R  -  CP  *  °V 


Y  =  cp/cv 


(7) 

(8) 

(9) 

(10) 


Hence  we  may  put  Equations  (1),  (2),  and  (3)  into  the  form 


—  -  —  (1  -  §) 

p,  p  0 

I  o 


P,  =  P  +  p  DU 
f  o  o 


p  ,  P,  +  P  i  . 

O  f  f  O  ,  1  1  V 

p  (y-D  +  P  o  +  2  >  "  p  ; 

o  o  or 


pf(y-D 


(ID 

(12) 

(13) 


The  pressure  density  states  achievable  for  a  given  f  are  those  given  by 
Equation  (13)  where  D  is  to  be  eliminated  through  (11)  and  (12).  Combining 
(11)  and  (12)  to  eliminate  U  yields  (using  n  =  PQ/Pf  as  8  compression  factor) 


U  =  D(l-n) 

Pf  -  P  +  p  D*(l-n) 

too 


(14) 

(15) 


and  finally 


*  (pf  *  po)(n 


-  1) 


(16) 


PQ(l-n) 
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This  is  readily  solved  for  if  we  assume 

PQ  «  Pf  (17) 

whence 

pf3/2  =  2f  ✓po(l-n)  "  l]  <l8) 

This  result  has  been  given  by  Razier' who  argues  that  the  steady  detonation  wave 
occurs  when  the  sonic  speed  behind  the  front  given  by 


c 


t 


(19) 


is  equal  to  the  shock  propagation  speed,  measured  with  respect  to  the  wave. 
The  sonic  speed  from  (15)  and  the  definition  of  n  yields 


c*  =  YD*(l-n)n 


Thus,  the  condition 


c  =  D  -  U 


yields,  after  using  (14),  that 


Taking  this  to  Equation  (18)  yields 


(20) 


(21) 


(22) 


[ 


2f  (y-1) 


(23) 


^Raster,  Y.  P.  "Heating  of  a  Gas  by  a  Powerful  Light  Pulse,"  Soviet 
Physics  JETP.  Vol.  21,  Ho.  5,  1965,  pp.  1009-1017. 
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Substituting  (22)  into  (15)  gives  the  detonation  speed  in  terms  of  pressure  as 


i  D 

f  y+1 


from  which  we  obtain  Razier's  result 


D 


(y2-D 


(24) 


(25) 


for  the  detonation  speed  as  a  function  of  energy  absorbed. 

Thus,  we  expect  a  cylinder  of  air  to  be  shocked  to  this  pressure  P^,  the 
length  of  the  cylinder  increasing  with  time  as  the  detonation  front  propagates 
up  the  beam  at  speed  D. 

However,  the  high  pressure  in  the  cylinder  of  air  will  relax  as  the 
cylinder  expands  radially  as  a  cylindrical  blast  wave.  Initially  the  cylinder 
is  of  beam  radius  a,  and  pressure  P^,  given  by  (23).  As  the  cylindrical  blast 
wave  expands,  the  same  Rankine  Hugoniot  Jump  relations  must  be  satisfied  across 
the  shock,  which  we  take  to  be  strong,  i.e. 


(26) 


Hence  the  radial  particle  velocity  is  related  to  the  shock  speed  by 


U 

r 


=  R(1  - 


Y+r 


R 

Y+l 


and 


(27) 


(28) 
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From  considerations  of  similarity,  as  in  the  Taylor  Blast  Wave  solution, 
we  note  the  energy  per  unit  length  in  the  s  direction  of  the  expanding 


eye  Under  is 


rR  P(C)P  R*  •, 

*  fitly— 1^1  *  T  U‘<?»  2,rdr 

o  t 


where 


Hence 


must  be  constant,  or 


from  which  we  deduce  that 


5  *  r/R 


P(r)  =  pqR*PU) 


U(r)  *  RU(C) 


■  \  ({&  *  1 2,R‘S' 


RR  =  2C* 


R  =  2C  /t  |  R  =  C(/t)_1 


must  describe  the  growth  of  the  cylindrical  blast  wave. 

For  an  Initial  radius,  a,  and  initial  pressure,  P^,  we  see 


Hence 


(Y+DP, 


/ (?+l)Pf 

2c*  =  a/  __I 


=  a(y+l)c. 


where 


Finally, 


IHP 

i  °t 


R*R*  =  a*Pf(1±i) 
f  po 


Pf 

P(R)  =  a1 


R  =  <2a(y+l)cft)i/2 


This  result  differs  slightly  from  that  of  Nielsen  who  writes 


R  ■  /2c  ak 
0  t 


Nielsen,  P.  E.  "Hydrodynamic  Calculations  of  Surface  Response  in  the 
Presence  of  Laser-Supported  Detonation  Waves,"  Journal  of  Applied  Physics. 
Vol.  46,  No.  10,  October  1975,  pp.  4501-4505. 
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where 


From  the  foregoing,  we  may  expect  that  a  high  intensity  laser  will  produce 

a  pressure  pulse  on  a  surface  placed  in  the  beam.  The  pressure  will  decay 

quickly  as  a  blast  wave  spreads  radially.  For  high  Intensity  lasers  operating 

in  a  pulsed  mode,  the  result  will  be  a  pressure  field,  varying  both  spatially 

and  temporally  in  a  very  complex  manner. 

The  mechanism  of  pressure  generation  we  have  considered  is  that  known  as 

the  LSD  (laser-supported  detonation)  wave.  At  lower  intensities,  a  different 

3 

process,  the  LSC  (laser-supported  combustion)  wave  occurs  .  In  this  case,  the 
energy  balance  on  the  processed  column  of  air  (plasma)  must  include  the 
radiation  loss  to  the  target  material.  Numerical  studies  have  shown  the 
pressure  generated  by  this  mechanism  to  be  50-60  atm,  with  up  to  50X  of  the 
incident  flux  being  re-radiated  to  the  target.  The  absorption  of  this  radia¬ 
tion  produces  the  process  known  as  thermal  enhancement. 


No. 


^Thomas,  P.  D.  "Laser  Absorption  Wave  Formation,"  AIAA  Journal, 
10,  October  1975,  pp.  1279-1286. 


Vol . 


13, 
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B.  The  Response  of  Finite  Plates 

The  response  of  an  elastic  plate  to  a  spatially  and  time  varying  distri¬ 
bution  of  pressure  arises  in  several  contexts  of  the  general  problem  of  deter¬ 
mining  the  response  of  structures  to  High  Intensity  Laser  radiation.  The 
differential  equation  of  the  elementary  theory  of  plate  bending 

dv*w  +  ph  |pr  =  P(*>y>t)  O) 

where  w  is  the  displacement  normal  to  the  plane  was  encountered  earlier  in 
determining  the  response  to  a  thermally  induced  distributed  moment.  Here,  we 
restrict  our  attention  to  the  isothermal  deformations  of  a  simply  supported 
circular  plate  under  axisymmetric  loading.  Equation  (1)  is  obtained  by  noting 
that  Equation  (9)  (of  Section  II-C)  must  be  replaced  by 

fr  -  fr 

2wrQ  =  j  ph2irw(r,t)rdr  -  I  P(r,t)2irrdr 

o  o 

where  P(r,t)  is  an  axisymmetric,  distributed  pressure.  Equations  (5a)  and  (5b) 
go  through  with  the  only  change  being  the  deletion  of  the  thermal  moment.  For 
the  simoly  supported  plate,  the  eigenfunctions  are 

R(r)  =  J  (pr)  “  J  (p  a)I  (p  r)/I  (p  a)  (2) 

n  o  n  o  n  o  n  o  n 

with  eigenvalues  satisfying  Equation  (II-E-24).  The  first  several  are  given 

1  2 

(for  v  =  .3)  in  Leissa's  monograph  and  others  in  more  recent  work  .  The  first 


^Leissa,  A.  W. ,  Vibration  of  Plates,  NASA  SP-160,  1969. 

2 

Lcissa,  A,  W.  and  V.  Narita,  "Natural  Frequencies  of  Simply  Supported  Circular 
Plates,"  accepted  for  publication.  Journal  of  Sound  and  Vibration. 
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11  ara  tabulated  In  Table  1  for  2  values  of  Poissons  ratio.  As  each  eigenfunc¬ 
tion  satisfies  the  condition  of  vanishing  moment  and  displacement  on  r  =  a, 
any  superposition  does  as  well,  and  the  general  solution  to  Equation  (1)  may 
be  written  as 


w(r,t)  =  l  W  (t)R  (r) 
,  n  n 


(3) 


We  expand  the  prescribed  pressure,  P(r,t)  as  a  series  of  these  eigenfunctions 


P(r,t)  =  l  f  (t)R  (r) 
n=l  n  " 


(4) 


Substitution  into  Equation  (1)  leaves,  after  multiplication  by  R  (rj  ar.a  inte- 

n 

gration  over  the  domain,  a  result  analogous  to  II-C-33,  or 


Dp'V  +  phVn  =  f  (t) 
n  n  n  n 


(5) 


where 


P(r,t)R  (r)rdr  . 
n 


Solutions  to  Equation  (5)  must  satisfy  desired  initial  conditions;  typically, 
that  the  plate  is  at  rest  for  negative  time  and  is  set  in  motion  by  a  pressure 
distribution  first  achieving  nonzero  values  for  t  =  0. 

As  a  first  example,  consider  a  suddenly  applied  uniform  pressure  F^,  i.e., 

P(r,t)  =  P  u(t)u(b— r)  (6) 

o 


Here  u(s)  is  the  unit  step  function 
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u(s)  =  0 

for  s  <  0 

u(s)  =  1 

for  s  >  0 

(7a) 

(7b) 


Then  ^n(t)  =  0  for  t  <  0,  while  for  all  t  >  0 


£n(t)  *  £„°  '  5®  f  Vr)rdt  *  f 

n  o 


K  J  (P  »)  h 

-  J.(p  b)  -  - — - r  —  I.(p  b) 

<>n  1  "  V».*>  Pn  1  “ 


(8) 


and 


=  a* 


fV 


-  J  *<7^ 
o  1-v 


-  I  */I  *) 
2  A1  '  o  ' 


] 


The  modal  amplitudes  azo. 


W 


n 


sin  w  t  +  B 
n  n 


COS  ti>  t  + 
n 


(9) 


(10) 


and  satisfaction  of  the  homogeneous  initial  condition  requires  that 


B 

n 


A  =  0 
n 


(ID 


Hence 


w(r,t)  = 


p  b  -  <Vpn»>>  -  riFTT  Il(pnb)} 

A  A  rfl 


-2-  l 

D  L. 
n=l 


o  n 
[Q(*  V 
Nn  Pn 


{Jo(pnr) 


J  (p  a) 

T°/  ~  t  1  (p  r'JHl  -  cos  (i)  t)  (12) 
I  (p  a)  o  .»  n 


where 


(13) 
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Fn(b/a,n) 


J  (pa)  J  (p  a) 

ui<»„b>  -  tVjt  Ii<p»b),<1  * 

_ o  rn  o  n  _ 

.  x(p  a)  ,  ,  l/(p  a) 

—  2 - Jo  (pna)  <U^  +  2  I  r(p  a)}  (pna) 

o  n 


M  =  P  af(b/a)(— tt)  l  F  (b/a,n)(l  -  cos  u>  t) 
max  o  z  .  n  n 

n=i 


and  the  maximum  stresses  are 


H 

a  =  3(l+v)P  (—)*(“)  l  F  (b/a >n)  ( 1  -  cos  u  t) 
max  o  n  a  .  n  n 

n=l 


Table  II  presents  values  of  for  several  b/a  and  several  modes  for  material? 
with  v  =  .3,  and  may  ba  used  to  determine  which  modes  contribute  most  signifi¬ 
cantly  to  the  stress  history  for  any  b/a.  If  b  =  a,  the  first  mode  predominates 
and  we  find  (for  v  =  .3) 


F^l.l)  =  .344 


The  stress  history  at  the  plate  center  is 


a  =  1.34  P  (a/h  )  *(1  -  cos  w.t) 
max  o  1 


1.  Response  to  a  Single  Pulse 


We  may  deduce  the  response  to  a  square  pulse, 


P  =  P  0  <  t  <  T 
o  -  — 


(24a) 


P  =  0  t  >  T 


(24b) 
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p 


TABLE  II  F  (b/a,n)  for  v  =  .3 

VI 


b/a=. 1 

b/a=.3 

b/a=.5 

.07601 

.21608 

.32204 

.02802 

,06120 

.04644 

.01666 

.02021 

-.00739 

,01120 

.00286 

-.00679 

.00788 

-.00316 

.00175 

.00561 

-.00333 

.00234 

.00396 

-.00135 

-.00070 

.00272 

.00044 

-.00111 

.00178 

.00109 

.00036 

.00106 

.00074 

.00063 

.00052 

.00003 

-.00021 

TABLE  III  F  (b/a,n)  for  v  -  .2* 
n 


b/a=. 1 

b/a=.3 

b/a= .  5 

.07903 

.22482 

.33556 

.02815 

.06156 

.04693 

.01669 

.02020 

-.00737 

.01121 

.00287 

-.00681 

.00788 

-.00316 

.00174 

for  n  ^  6,  values  same  as  for  v  =  .3 


b/a=l  .0 

.3*395 

-.03441 

.01113 

-.00515 

.00286 

-.00178 

.00119 

-.00085 

.00063 

-.00048 

,00037 


b/a=l  .0 

.35990 

-.3431 

.01109 

-.00513 

.00285 


by  superposition.  For  t  <  i ,  Equation  (21)  applies.  For  t  >  t,  we  find 


o  =*  3(l+u)P  (r)*(~)  I  F  (b/a,n)[cos  m  (t-t)  -  cos  w  t] 
on  a  n  n  n 


The  maximum  contribution  for  any  one  mode  occurs  when  t  =  t^,  where 


tan  a)  t  =  -cot(u>  t/2) 
n  c  n 


A  useful  approximation  results  for  pulses  of  duration  much  less  than  the  period 


of  the  dominant  modef 


u  t  =  n/2 
n  c 


where  the  impulse  per  unit  area,  Iq  =  Pqt,  and  N  denotes  the  dominant  mode. 
For  a  uniform  pressure,  the  dominant  mode  is  the  first  and 


o  =  1.34  I  (a/h)*w. 
max  o  i 


For  a  nonuniform  pressure,  the  first  mode  response  gives 


V  “l^ot.l  ,  F1  > 
— hi —  lb7I> 


Recall  from  Equation  (13)  that 


u  =  J  ”T  (p  a)1 
n  v  Ph  a*  *n 


TI(lW  (£)(pna)‘ 
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Thus,  Che  necessary  modes  for  the  modal  analysis  are  governed  by 

(p  a)*F  (b/a,n). 
n  n 

This  result  may  be  used  for  comparison  with  computations  performed 
by  Holmes,  Keough  and  Desmond^  for  pyrocerams  (E  =  16.5  x  10*  psi; 
p  =  .0941  lb/in*)  v  =  .25,  a  =  4  in.,  h  =  1/4  in.  subjected  to  a  pulse 


P(r,t>  .  P 

0 

t  =  20  usee,  s  =  2.75  inches,  Pq  =  1087  psi 

Such  a  pressure  distribution  history  gives  a  total  impulse  of 

(32) 

Itotal  -  f  2»(P(r,t)  )rdr  dt 
o 

(33a) 

-  tP  2tt  f  e-(r/s),rdr 
°  Jo 

(33b) 

=  TP  ws*(l  -  e-U/s)‘) 
o 

(33c) 

=  (2.17  x  10-2  psi  sec)xs* ( .8795) 

(33d) 

which  is  equivalent  to  a  uniform  pressure  Pq  over  b  =  .938  s. 

For  this  disk, 

a)  =  6.127  x  10*  rad/sec 

(34) 

and  the  period  is 


T 


1.025  x  10 


sec  =  1025  usee 


(35) 


We  obtain  a  maximum  stress  from  Equation  (30) 


^Holmes,  B.  S.,  D.  D.  Keough  and  T.  P.  Desmond,  Monthly  Status  Report  #2, 
SRI  International  Project  PYU-7259,  Contract  No.  F29601-78-C-0041 ,  August  1978. 
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=  3.9(2.17  x  10"2  psi  sec)(. 8795X6.127  x  10*/sec) 

.  F  (.645,1) 

•  <“t)  ■  -  ;643 -  =  30.6  ksi  (36) 


occurring  at  the  quarter  period  of  the  first  mode,  256  usee. 

Values  of  the  function  Fn(b/a,n)  and  the  product  (pna)*Fn(b/a,n)  are  given 
in  Tables  III  and  IV  for  v  =  .2.  Neither  is  very  sensitive  to  the  value  of  v. 
The  function  Fr  was  seen  from  Equation  (21)  to  be  useful  in  determining  the 
dominant  mode  in  the  response  due  to  the  step  application  of  a  uniform  pressure, 

a 

and  the  tabulated  values  show  that  the  dominant  mode  is  the  first.  The  values 
tabulated  as  Table  IV  demonstrate  that  the  response  to  a  pulse  loading  which 
does  not  flood  the  entire  plate  will  contain  strong  contributions  from  a  number 
of  modes.  For  example,  in  the  case  of  b/a  =  .1,  the  contribution  of  the  11th 
mode  is  greater  than  the  first. 


2.  Response  to  a  Train  of  Pulses 

Let  us  consider  the  response  of  a  simply  supported  elastic  plate  of  radius 
a  to  a  train  of  pulses  of  a  particular  temporal  variation: 

P(r,t)  =  Pq  U(b-r)[l  -  cos  tit]  for  t  >  0  (37a) 

P(r,t)  =0  for  t  <  0  (37b) 


This  may  be  considered  to  approximate  a  train  of  "square"  pulses,  of  duration. 
T/2  and  amplitude  2Pq,  recurring  with  period  T  =  2*/n,  as  shown  in  Figure  9. 
Each  pulse  delivers  an  impulse 


Ip  —  |  |  P(r,t)2itrdr 

1  o  ‘  o 

#2*/n 

I_  *  iP  b*  (1  -  cos  Ot)dt  =  P  wb*T 

r  O  ]  o 

o 
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(38a) 


(38b) 


i  i  t"  ■  V  i<l*Yi*ra<)  i  i  *  i-  r  i 


2 


TABLE  IV  (p  a)*F  (b/a)  for  v  =  . 
n  n 


b/a=.l 

b/a=.2 

cr 

ll 

• 

to 

b/a=.5 

b/a=1.0 

n 

I 

.37797 

.74113 

1.07522 

1.60483 

1.72125 

2 

.83325 

1.48629 

1.82233 

1.38935 

-1.01556 

3 

1.23567 

1.82646 

1.50143 

-  .54562 

.82101 

4 

1.54893 

1.66943 

.39704 

-  .94082 

-  .70896 

5 

1.75084 

1.08918 

-  .70196 

.38713 

.63343 

6 

1.82879 

.28746 

-1.08477 

.76221 

-  .57789 

7 

1.77999 

-  .47691 

-  .60673 

-  .31562 

.53485 

8 

1.61179 

-  .97354 

.26017 

-  .65797 

-  .50023 

9 

1.34082 

-1.07376 

.82247 

.27297 

.47162 

10 

.99163 

-  .78715 

.69000 

.58750 

-  .44732 

11 

.59410 

-  .24977 

.03482 

-  .24383 

.42650 

69 


ia aea 


► 

Li 


The  modal  amplitudes  are  given  by 
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As  beforei  we  find  from  Equation  (9),  and  that 
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We  first  note  the  possibility  of  resonance  if 


It  is  to  be  noted  here  that  in  the  pulse  train  modeled,  T  is  the  pulse  spacing, 
rather  than  the  pulse  duration. 

For  a  long  train  of  pulses,  with  T  near  to  T  , 
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(3) 

For  long  times,  Wn  will  cease  to  be  important,  through  energy  dissipation  not 
included  in  this  analysis.  Ve  have,  then,  for  such  long  times 
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where  the  eigenfunctions,  Rn(r),  are  as  given  previously.  Equation  (2).  We  may 
again  compute  moments,  ond  stresses,  finding  the  maximum  stress  at  the  plate 
center  (for  long  times)  to  be 
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At  short  times,  a  greater  stress  may  occur  at  t  =  such  that 

(n*/<o  1  -  2)  sin  w  t,  =  —  sin  fit, 
n  n  1  u>  1 

n 


(53) 


of  magnitude 


max 
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If  II  **  the  forced  response  [Equation  (54)]  may  become  very  large.  If,  in 
fact,  $1  =  d»N ,  we  find 
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It  is  not  likely  that  a  target  can  be  driven  precisely  at  one  of  its  resonant 
frequencies;  however,  if  this  could  be  achieved,  the  stress  level  after  M  pulses 
would  be  (for  a/b  =  1  and  resonating  the  first  mode) 
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For  the  pyroceraa  disk  previously  cited,  a  stress  of  18  ksi  would  be  achieved 

_2 

in  100  pulses  of  .7  x  10  lb-sec  per  pulse,  as  opposed  to  the  single  pulse  of 
about  60  times  that  magnitude.  That  the  resonance  effect  is  not  more  effective 
is  somewhat  surprising.  The  explanation  probably  lies  in  the  fact  that  the  long 
train  of  small  pulses  produces,  even  at  resonance,  what  is  essentially  a  static 
response,  whereas  the  single  pulse  of  short  duration  allows  the  dynamic  ampli¬ 
fication  (typically  a  factor  of  two)  that  is  characteristic  of  impulse  response. 
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IV.  A  Fracture  Criterion  for  Brittle  Materials 


A.  Introduction 

Experiments  have  shown  that  a  ceramic  material  heated  over  some  portion 
of  its  surface  will  fracture.  These  fractures  appear  to  originate  at  the  rear 
surface,  under  the  heated  zone.  Consequently,  we  deduce  that  radial  conduc¬ 
tion,  even  in  poor  conductors  such  as  ceramics,  is  sufficient  to  prevent  large 
tensile  stresses  from  forming  outside  of  the  heated  region,  as  described  in 
Section  II. 

t 

Two  other  investigators  have  developed  criteria  for  the  fracture  of 
ceramics.  Mechlosky^  assumes  that  stresses  of  magnitude 


«  -  E  <*V  U) 

are  developed  and  that  fracture  occurs  when  a  critical  value  is  reached.  He 

2 

takes  the  temperature  rise  to  be  that  of  Goodmans  premelting  approximation 


where 


t' 


(3) 


Solving  for  the  time  t  to  reach  a  critical  stress,  o  ,  we  find 

c 


tet  1  1  /,  xs,  2k  (l-v)o 

h*  =  e  "2 (1  -h0  +7T^r~ 

a 


(4) 


^Mecholsky,  T.  T.,  et  al.,  "Effect  of  3.8u  Laser  Radiation  on  Ceramics", 
Proceedings  of  2nd  Conference  on  Laser  Effects  and  Hardening,  NASA  Ames, 

July  1975. 
o 

Goodman,  T.  R.  and  J.  J.  Shea,  "The  Melting  of  Finite  Slabs",  Journal 
of  Applied  Mechanics,  Vol.  27,  pp  16-24,  March  i960. 
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A  factor  of  two  has  been  introduced  into  the  last  term  to  account  for  the 
radial  expansion  of  the  surrounding  material,  as  in  Section  II,  and  is  needed 
to  produce  Mechlosky's  result. 

Mechlosky  assumes  that  the  critical  value  of  x  is  h/2,  thus  he  finds  a 
fracture  time 


Ktf  _  _L  _2k_  (l-v)o  ,*) 

P"  “  24  +  I  h  cxE  V 

a 

For  sufficiently  high  flux  (low  strength)  the  right  hand  term  becomes 
negligible  leading  to  a  dimensionless  fracture  time  which  is  Independent  of 
intensity. 

Several  observations  regarding  this  criterion  must  be  made.  First,  the 
stress  distribution  assumed  (Equation  1)  is  appropriate  only  for  predicting 
the  compressive  stress  arising  under  two-dimensional  constraint,  or,  with  the 
inclusion  of  the  factor  of  2,  partial  constraint.  No  out  of  plane  deformation 
is  allowed.  Stresses,  in  this  model,  are  compressive  everywhere  and,  as  such, 
the  stresses  at  the  midplane  would  appear  to  have  no  particular  significance. 
Further,  the  model  does  not  account  for  the  observed  origin  of  fracture  at  the 
rear  surface.  Finally,  the  result  is  extremely  sensitive  to  the  value  of  x 
employed.  Nonetheless,  it  must  be  noted  that  quite  good  results  have  been 
obtained  with  it. 

3 

Laughlin  has  developed  a  means  of  estimating  fracture  time  based  on  the 
assumption  that  fracture  will  occur  when  tensile  stresses  on  the  rear  face 
reach  a  critical  va.uc.  He  begins  by  using  the  exact  solution  for  the 


3 

Laughlin,  W. ,  "Predicting  the  Laser  Induced  Thermal  Fracture  of  Infrared 
and  Radar  Transmitting  Materials",  Proceedings  of  First  Classified  Conference 
on  High  Energy  Laser  Technology,  San  Diego  CA,  October  1974. 
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temperature  rise  In  a  uniformly  heated  slab  to  predict  thermal  stresses, 

under  the  assumption  that  the  material  surrounding  the  heated  area  constrains 

all  bending  deformation  but  allows  radial  expansion  to  freely  occur.  This 

corresponds  to  the  mechanical  boundary  condition  discussed  as  case  II  in 
'  4 

Section  II. B. 2.  To  account  for  melting  ,  he  computes  stresses  in  the  same 
manner,  with  the  temperature  distribution  taken  to  be  that  in  a  slab  with 
front  face  temperature  fixed  at  meeting.  As  the  slab  continues  to  heat, 
we  note  that  the  front  face  gradient  decreases,  and  thus,  the  incident  flux 
is  treated  as  if  it  were  reduced.  Laughlin  has  also  obtained  generally  sat¬ 
isfactory  results  with  his  method. 

While  Laughlin's  model  appears  to  have  a  firmer  basis  that  that  of 
Mechlosky,  two  questions  remain  open.  First,  the  mechanical  boundary  con¬ 
dition  at  the  edge  of  the  heated  zone,  and  secondly,  the  treatment  of  the 
influence  of  melting  on  the  stresses  in  the  solid  portion. 

B.  The  Influence  of  Mechanical  Boundary  Condition 

Let  us  consider  a  section  of  a  structural  member  to  be  subjected  to  a 
uniform  absorbed  intensity  over  a  circle  of  radius  a,  the  remainder  of  the 
structure  remaining  at  the  ambient  temperature  which,  for  convenience,  will 
be  taken  as  T  =  0.  If  we  neglect  the  consequences  of  radial  heat  flow  out 
of  the  heated  region  ,  then  the  temperatures  for  r  <  a  are  given  by  Equa¬ 
tion  (18)  of  Section  II. B. 2.  or,  in  terms  of  a  coordinate  system  in  the 
midplane  of  the  plate  with  the  heated  face  being  z  =  h/2,  by 


Recall  that  it  was  shown  in  an  earlier  section  that  radial  heat  flow 
does  not  reduce  the  magnitude  of  stresses  within  the  heated  zone. 

4 

Laughlin,  W.  T.,  "Predicting  the  Laser  Induced  Thermal  Fracture  of 
Partially • Transparent  Materials  (U)",  2nd  L'OD  Conference  on  Laser  Effects 
and  Hardening,  NASA  Ames,  July  1975. 
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(1) 
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In  what  follows,  a  dimensionless  time 

t  =  iet/h2 


(2) 


will  be  used. 

For  r  <  a,  the  displacements  and  stresses  are  found  to  be 

o  =0 
zz 

o  —  o  =  i"  [C,Z  +  C*  —  aT  J 
rr  86  1-v  L  1  2  J 

U  =  C  Zr  +  C„r 
r  1  2 

U0  =0 

where  the  constants  and  are  to  be  determined  from  the  boundary  condi¬ 
tions  on  r  =  a.  The  four  possible  boundary  conditions  arising  from  the  edge 
r  =  a  being  free  or  constrained  against  in  displacement  and  free  or  con¬ 
strained  against  rotation  were  explored  in  a  previous  section. 

We  now  assume  the  line  r  =  a,  circumscribing  the  heated  zone,  to  be 
elastically  restrained  by  the  surrounding  material,  i.e.,  on  r  =  a  the  in¬ 
fluence  of  the  surrounding  material  may  be  represented  by  moments  and  forces 
per  unit  length 
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Further,  the  displacements  and  rotations  of  the  heated  region,  on  the  line 
r  =  a,  must  match  those  of  the  exterior  region.  Here,  Ur  and  w  are  the 
inplane  and  out  of  plane  displacements,  respectively  in  the  exterior 
(unheated)  region. 

The  constants  of  (3b)  and  (4a)  result  from  matching  displacement,  slope, 
force  and  moment  at  r  =  a 
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(11c) 
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These  four  equations,  together  with  the  two  stiffness  constants,  defined  by 
Equations  (5)  and  (8),  permit  the  constants  and  to  be  determined.  We 
find: 
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fh/2 

a  Tdz 

-h/2 
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E  T 


1  fh/2 
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(12) 


(13) 


The  four  cases  considered  in  Section  II. B. 2.  may,  of  course,  be  recovered  by 
setting  to  zero,  or  infinity,  the  two  stiffnesses,  CT  and  C ,  as  appropriate. 

We  first  consider  the  flat  plate  of  thickness  h  surrounding  the  heated  re 
gion.  We  provide  that  the  elastic  constants,  E,  v,  may  differ  from  those  of 
the  heated  region  so  that  a  thermal  softening  effect  may  be  included,  if 
desired.  The  inplane  displacements  and  stresses  are  those  of  the  Lame  solu¬ 
tion,  as  used  previously  (Section  II-B). 
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(15b) 
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For  a  finite  circular  plate,  with  Ff  =  0  on  the  outer  radius,  r  =  b,  we  find 


B  .  A(|±J)b 


(16) 


For  an  infinite  sheet,  A  =  0. 

These  results  may  be  combined  to  deduce  the  stiffness  constant  for  the 
inplane  motion,  C^.  Using  Equations  (8),  (9),  and  (10),  we  find 
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In  order  to  facilitate  comparison  with  a  fracture  criterion  developed  by 
Laughlln,  it  is  of  interest  to  determine  the  stresses  arising  in  a  finite 
sheet  of  the  same  modulus,  with  the  central  region  completely  constrained 
against  bending,  Cg  =  »  .  In  this  case, 
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From  Equations  (12)  and  (13),  =  0  and 
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Since  the  moment  et  r  *  b  oust  vanish,  the  moment  at  r  =  a  is  found  to  be 
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The  slope  at  i  =  a  is  seen  to  be 
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from  which  Equation  (7)  may  be  used  to  deduce  the  stiffness  as 


c  2h>  (1-^ 

B  12(l+v)a  7i  1 -v  a 1 „ 

(1  +  it:  bT) 


Hence 


rh/2 

zTdz 

-h/2 


[1+f  (1-v) 


(1+v  +  (1-v)  p-) 


In  the  case  of  an  infinite  sheet  with  elastic  properties  on  r  <  a,  these 


reduce  to 
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In  order  to  evaluate  the  stresses  from  Equation  (3b),  we  need  compute 


certain  integrals  of  the  temperature  distribution 
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Let  us  examine  the  following  cas.es: 

A:  The  edge  r  =  a  is  completely  constrained  against  rotation,  but  is  only 
partially  constrained  against  lateral  expansion  by  the  radial  deformation  of 
a  finite,  circular  plate.  The  constant  C^,  as  noted  previously,  is  zero  and 

(1+v)  +  (1-v)  fh/2 

c  =  [ - - — ]  t-  Tdz  (35) 

h  '-h/2 


The  stresses  in  r  <  a  are 
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For  this  boundary  condition,  the  maximum  stress  occurs  on  the  rear  face 
(*  =  -h/2)  and  is 
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(37) 


Figure  1  displays  graphically  this  equation  for  three  cases:  a/b  =  1  (which 


corresponds  to  Laughlins  criterion);  a/b  =  0  (which  corresponds  to  the  infi¬ 
nite  sheet  with  a  heated  spot  and  an  intermediate  case  (a/b  =  1/2),  all  with 
v  =  .3. 

Qualitative  differences  between  the  two  extreme  cases  are  seen.  For 
a/b  =  1,  the  rear  face  stress  increases  monotonically  reaching  a  maximum 
(as  t  *  •)  of 


En  Ia*>  1,  Ea  *a** 
°max  -  1-v  k  X  6' 1-v  k 


(38) 


On  the  other  hand,  the  elastic  constraint  caused  by  the  infinite  sheet  leads 
to  the  eventual  (t  -  1/2)  development  of  compressive  stresses  on  the  rear 
face,  with  the  maximum  tensile  stress  occurring  at  t  =  .175,  of  magnitude 
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From  the  results  given  in  Figure  IQ  a  failure  time  can  be  predicted.  If  the 
case  a/b  =  1  is  used,  Laughlins  rer «its  will  be  recovered.  It  is  interesting 
to  note  that  for  times  in  the  region  where  most  failures  have  been  observed, 
that  the  case  for  a/b  =  l  overestimates  the  stress  by  about  50%  or  equiva¬ 
lently  provides  an  estimate  of  the  fracture  time  which  is  about  60%  of  that 
which  would  be  obtained  if  the  a/b  =  «•  results  were  to  be  used. 

Let  us  now  consider  a  second  case:  where  the  line  r  =  a  is  only  par¬ 
tially  constrained  against  both  expansion  and  rotation  by  the  surrounding 
material.  C2  is  as  before,  but  is  no  longer  zero,  for  the  bending 
deformation  will  now  take  place.  Using  (33)  and  (34)  in  (35)  and  (30) 

(taking  E  =  I  and  v  =  v)  we  find  the  stress  (Equation  (3b)  to  be 
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where 
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For  a  =  b,  R  =  2,  and  we  find,  for  long  tines, 
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which  is  the  result  previously  obtained  for  the  free  plate.  For  b  =  ®,  the 
long  time  stresses  are 
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Evaluating  Equation  (40)  at  the  rear  face  (s  =  -h/2),  we  find  that 
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which  is  seen  to  differ  from  Equation  (37)  (the  result  when  bending  deforma 
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Stresses  at  the  rear  face,  for  short  times,  are  given  in  Figure  11  and  are  seen 
to  remain  compressive.  Long  time  values  (Equation  (42)  and  Equation  (43))  are 
also  seen  to  be  compressive. 

Thus,  we  see  that  properly  accounting  for  the  mechanical  boundary  con¬ 
dition  at  the  edge  of  the  heated  region  shows  the  stresses  on  the  rear  face 
to  remain  compressive. 

Tensile  stresses  will  develop  in  the  interior.  It  is  to  be  recalled  from 
Chapter  II  (Figure  6)  that  a  uniform  flux  over  one  face  of  an  infinite  sheet 
leads  to  compressive  stresses  on  the  outer  surface  and  tensile  stresses  on 
the  centerline.  These  results  are  reviewed  as  Figure  12a.  if  we  now  consider 
a  heated  radius  a  within  a  disk  of  outer  radius  b,  the  constraint  of  the 
unheated  area  superposes  compressive  stresses,  leading  to  the  time  history 
shown  in  Figure  12b  (for  a/b  -*•  0).  Small  compressive  stresses  remain  on  the 
center  only  for  short  times. 

Thus,  if  fracture  due  to  tensile  stresses  on  the  rear  surface  is  observed, 
the  fracture  must  originate  outside  the  heated  area,  according  to  the  mechanism 
of  Chapter  IIB. 

The  stresses  in  a  disk  of  radius  b  having  an  axisymmetric  temperature 
distribution  are  given  in  Equations  II-A-19  and  20.  For  a  temperature  distri- 
but ion 


T  = 

T 

0 

on  r  _<  a 

(46a) 

T  = 

0 

on  a  <  r  <  b 

('‘6b) 
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Figure  12.  Stress  Distribution  at  Several  Times 


we  find 


o..  =  a  =  EaT  (£t  -  l)/2  on  0  <  r  <  a  (4?a) 

88  rr  o  D 

a  =  EaT  (rr  ~  t*)/2  on  a  <  r  <  b  (47b) 

rr  o  b1  r*  —  — 


and 


-  E“To(b 


*  «* 

+  7t)/2 


on 


<  r  <  b 


(48) 


Thus,  the  greatest  tenslie  stress  occurs  at  r  =  a,  just  outside  the  heated  area. 

For  no  losses,  and  no  radial  conduction,  the  average  temperature  on  r  £  a 
is 


(49) 


and  the  corresponding  average  tensile  stress  for  r  -  a+  increases  linearly  with 
time  as 


a 

max 


Eal  2 
_ a  ,a 

=  TcT  (b*  + 

p 


(50) 


For  the  infinite  disk.  Equation  II-B-1-I2  is  recovered.  Radial  conduction  will 
lead  to  a  reduction  of  these  stresses,  in  the  manner  of  Section  II-A,  where 
asymptotic  value*-  for  long  times  have  already  been  given.  It  should  be  noted 
that  these  are  not  the  only  stresses  developed.  Continuity  of  bending,  as  well 
as  in-plane  deformation,  will  give  vise  to  a  moment  and  further  stresses 

o  .  These  will  not  be  considered  at  they  will  quickly  decay  with  distance 
into  the  unhested  portion.  Rather,  the  influence  of  radial  conductivity  on 
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tensile  stresses  for  r  >  a  will  be  explored.  The  stress  time  history  for  k  =  0 
Is  given  by  Equation  (50).  For  non-zero  tc ,  a  numerical  procedure  was  developed 
for  the  transient  on  dimensional  radial  temperature  distribution  in  a  finite 
disk.  Once  values  of  temperature  were  found,  stresses  were  computed  from 
Equations  II-A-19  and  20.  Typical  results  for  the  transient  circumferential 
stress  are  given  in  Figure  13.  Here  the  heated  radius  was  taken  to  be  4  cm  on 


a  disk  1  cm  thick  of  outer  radius  9.73  cm.  was  taken  to  be  unity  (1  Joule/ 
gra°K) ,  p  =  3  gm/cm*  and  two  values  of  conductivity,  leading  to  K  =  .1  cm* /sec 
and  SI  =  .01  cm* /sec,  were  considered.  The  radial  conduction  is  seen  to  reduce, 


but  not  eliminate,  the  development  of  the  tensile  stress. 


The  influence  of  conductivity  is  more  readily  seen  from  the  results  pre¬ 
sented  as  Figure  14.  A  family  of  hypothetical  materials  were  considered,  all 
with  pCp  =  3  Joule/ cm5 °K.  A  disk  9.73  cm  in  radius  and  1  cm  thick  was  con¬ 
sidered  to  be  subjected  to  a  constant  absorbed  flux  over  a  central  circle  of 


4  cm  in  radius.  The  time  to  generate  tensile  stresses  of  two  values  is 


presented.  From  the  results  it  can  be  seen  that  tensile  fracture  outside 


the  heated  region  may  indeed  occur,  and  that  the  time  to  produce  a  given  stress 
is  increased  by  thermal  conduction.  Thus,  heating  a  uniform  disk  will  produce 
significant  tensile  stresses  cutside  the  heated  region,  and  these  will  provide 
the  failure  mechanism  for  a  brittle  material. 
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Figure  13.  Tensile  Seres 


01  cmVsec 
1  cm2/sec 


Figure  14.  Fracture  Times  for  Hypothetical  Materials  of  Various  Conductivity 


C.  Influence  of  Melting 

When  the  heated  surface  of  the  slab  reaches  the  melting  temperature,  the 
temperature  distribution  is  no  longer  that  given  by  Equation  (II.B.1.),  and 
the  stresses,  as  computed  in  the  previous  section,  are  no  longer  applicable. 
Since  the  problem  of  determining  the  temperature  distribution  is,  in  such 
cases,  nonlinear,  we  must  resort  to  approximate  methods.  The  heat-balance 
Integral  method  of  Goodman  and  Shea  is  applicable  in  such  cases  and  will  be 
used  to  determine  an  approximate  solution  to  the  one-dimensional  heat  transfer 
problem  with  change  in  phase,  and  the  thermal  stress  will  then  be  determined 
from  this  approximate  solution.  It  will  prove  useful  to  recognize  two  crit¬ 
ical  distinctions,  first,  whether  the  slab  is  thin  or  thick,  and  second, 
whether  the  melt  is  completely  retained  or  completely  removed.  A  slab  will 
be  regarded  as  thick  if  fron“.  face  melting  begins  before  the  rear  surface 
temperature  changes  and  thin  when  the  converse  occurs. 

1.  Melting  of  Thin  Slabs  With  Melt  Retention 


For  temperatures  below  melting,  we  assume  the  temperature  distribution 

in  a  slab  of  thickness  J.  absorbing  a  flux  I  on  x  =  0  to  be 

a 


T  :•  a(t)  +  b(t)x  +  c(t)x* 

for  x  <  6 

(la) 

H 

II 

O 

for  x  >  S 

> 

(lb) 

where  6(t)  is  a  thermal  penetration  distance.  On  ::  =  0, 


on  x  -  5 


-kb(t) 


=  0  and  ~  = 
3x 


(2) 


(3) 
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This  approximation  cannot  satisfy  the  differential  equation  exactly  but  only 
in  an  average  manner.  Thus,  we  require 


r6  3  *T  f*  .  3T 

k  U?  =  pcP  n 
o  o  r 


(5: 


.  H  I  k  il  I 

K  Jx  .  3x 

6  o 


(6; 


which  may  be  integrated  to  give 


«(t)  =  /6xt 


(7 


since  S  =  0  at  t  =  0.  The  time  to  front  surface  melting  is  then  found  from 
(A)  to  be 


2kT 


(8 


2kT 

and  is  valid  if  6(t  )  =  — ; —  <  X.*  This  notation  suggests  the  introduction  of 
m  X 

a  dimensionless  measure  of  the  applied  flux.  Let 


i 


U 

~  (2kT  ) 
ro 


(9 


then  for  i  >  1,  melting  begins  before  the  rear  surface  temperature  changes, 
and  the  Blab  is  thermally  thick. 

We  first  consider  the  case  i  <  1,  the  case  considered  by  Goodman  and 


95 


2  l 

Shea  ,  with  results  used  by  Mecholsky  .  At  some  time  t.  less  than  t  ,  given 

•c  m 

by  (8),  the  penetration  depth  6  reaches  x  =  JL,  The  assumed  temperature  dis¬ 
tribution,  Equation  (4),  Is  no  longer  appropriate.  The  new  boundary  condi¬ 
tions  become 


3T 

I  =  -k-r~  on  x  =  0 
a  3x 


(10a) 


0  =  on  x  =  J. 

3x 


(10b) 


From  (  7 ) , 


(U) 


thus,  the  approximate  solution  for  t  >  t^  becomes 


Taia(t)  -£x  +  2ihxl 


(12) 


Again,  we  satisfy  the  averaged  differential  equation,  this  time  obtained  by 
integrating  through  the  entire  slab  thickness 


3T  I  .  3T  I  3  r  „ 

3x  i  "*  3x  3t  J  P 

*c  o  o 


Tdt  =  pC  a-4 
P  P 


(13) 


The  resulting  differential  equation 


a  = 


PC 


(14) 


is  easily  solved  and  substituted  into  (12)  to  yield  the  approximate  tempera¬ 
ture  distribution 
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(15) 


where  it  is  to  be  recalled  that  i  £  1. 

If  the  flux  is  continued  at  this  same  value  for  longer  periods  of  time, 
a  melting  front  will  form  and  propagate  into  the  solid,  leaving  behind  molten 
material,  with  thermal  properties  we  will  denote  by  a  subscript,  m. 

We  again  assume  parabolic  temperature  distributions,  measured  above  and 
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below  aelting,  1.#.,  T  ■  T  +  U  and  T  -  T  ♦  U,  respectively 


U  ■  d(t)  ♦  e(t)x  ♦  fft)x*  on  0  <  x  <  • 


U  «  a(t)  ♦  b(c)x  ♦  c(t)x*  on  a  <  x  <  A. 


(21a) 

(21b) 


with  boundary  condltlona  on  the  front  and  rear 


-  k 


3U 

a 

a  )x 


(22a) 


k|2  I  -0 
3X  U 


(22b) 


On  the  melt  Line,  x  =  s. 


U  =  U  =  0 

a 


(23) 


avid 


-  k 


3U 
_ m 

m  3x 


.  3U  I  ds 

*  k  v;  ■  PLdT 

s  s 


(24) 


The  first  four  equations  enable  us  to  reduce  each  of  (21)  to  two  unknown 
functions,  f(t)  and  C(t)  in 


U  =  —  (s  -  x)  +  f(t)(x*  -  s*) 
a  k 

ra 


(25) 


U  =  C(t)  (2/(s  -  x)  +  (x‘  -  s*)) 


(26) 


to  be  determined  by  the  heat  balance  Integral.  The  flux  matching  condition. 
Equation  (24),  then  provides  the  means  of  determining  the  unknown  time  hist^vy 
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of  the  panatratlon  depth,  s. 
Let 


e 

m 


•s 

U  dx 
m 
o 


s 


I 


r 


Udx  =  C(t)  l 


+  IU*  -  2aU 


|  s') 


Applying  the  heat  balance  Integral  on 


I 


s 

o 


3*U 


k 

m 


dx 


.s  3U 
_ m 

3t  m  pm 


dx 


and 


\l  k  0  d, .  rz  m 
L  3**  J„  3t 


P  C  dx 
in  pm 


produces 


i 


3U 

k  -r-S 

m  3x 


P  C 
m  pm 


de 

_ tn 

dt 


and 


-  k 


3U 

dx 


de 

dt 


pc 

p 


iut 


3U 

3x 


.  -  --o  - 
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and 


au 

IB 

ax 


.  I  38 

__L  +  f(t)29 

m  m 


(32) 


Substituting,  we  arrive  at 


_  de 
pCp  d£ 


3ke(t) 

“  TT^  ip 


(33) 


d8  31  3k  8 
_  in  ft  in  in 

Pra  pm  dt  ~  2  T* 


(34) 


*a  3kmVc)  3k8(t)  .  ds 

T  +  ~  Br~~  +  (7T if  =  fL  d7 


(35) 


At  the  onset  of  melting,  t  =  t  ,  s  =  0  and  6  =0.  The  temperature  distribu- 

in  n 

ticn  in  the  solid. 


U  =  C(t  ){x*  -  2ix) 
m 


(36) 


must  match  that  previously  given  by  the  prereelting  solution,  Equation  (16). 
Therefore 


C(t  ) 
in 


and 


8(t  )  =  - 
in 


i  i* 

a 

3k 


Introducing  the  same  new  variables  as  Goodman, 


(37) 


(38) 


(39) 
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(40) 


o> 


pc  e 

m  pm  m 

pC  T  l 
p  m 


v 


e 


=  T~L 

m 


T 


(t  -  t  ) 
m 


M  = 


C  T 
P  m 

L 


v 


k 

m 

< 


(41) 


(42) 


(43) 


(44) 


and  the  dimensionless  flux,  i,  as  previously  defined,  we  arrive  at 


dv  3v 

dt  “  “  (1  -  a)1 


(45) 


dw  , .  vu. 
dr  =  3(1 


(46) 


^  da  3va)  3y 

U  dx  +  a1  +  (1  -  a)* 


(47) 


with  o(0)  =  0,  ui(0)  =  0,  v(0)  =  -  —  these  are  the  same  equations  given  by 
Goodman,  who  has  given  approximate  solutions  obtained  by  a  perturbation 
technique. 

We  are  now  in  a  position  to  develop  expressions  for  Che  stresses  devel¬ 
oped  by  these  approximate  temper 'Cure  distributions.  Three  time  ranges  must 
be  considered.  To  simplify  comparison  with  Laughlin's  results,  we  will  con¬ 
sider  the  same  boundary  condition,  i.e.,  no  inplane  force  and  no  bending 
displacement.  Hence, 
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(48) 


o 


rr 


(C^x  -  aT) 


in  the  solid  region,  and  a  =  o.a  =  0  in  the  molten  region.  Hence, 

rr  00 


o 

rr 


ee 


Tdx  -  T] 


(49) 


Prior  to  the  time,  t^,  when  the  rear  surface  temperature  changes,  the  temper¬ 
ature  is 


with 


I  6 

T  =  — S—  (1 
2k  U 


for  0  <  x  <  6 


(50a) 


T  =  0 


for  6  <  x  <  J. 


(50b) 


6(t)  =  /6Kt  for  t  <  t  .  =  — 

-4  6k 


(51) 


a  --  c 

rr 


_  I  /6Kt  .  -TT7- 

Ea  a  i-l  /6Kt  ,,  x  ,  „  rrzr 

G00  -  Uv  2k  Cj  ~T  -  (l  -  ]  on  0  <  X  <  /6Kt 


(52) 


„  It 

Ea  a  .  . 

o  =  a  =  ■: - rr-  on  /oKt  <  x  <  J. 

rr  88  1-v  p^C 

P 


(53) 


Thus,  the  maximum  tensile  stress  is  at  the  rear  surface  and  rises  linearly  in 
time  reaching  a  maximum  value 


a 

max 


Ea 

1-v 


at  t 


(54) 


For  x 


>i 


but  before  the  front  surface  melts,  the  temperature  distribution 


is 
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(55) 
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is  reduced  to 


Eal  l 

o  =  [1  -  3(x  -  t  )]{T/na  ■;}  (61) 

max  *-  m  J  6k(l-v) 

Since  the  rear  surface  stress  depends  only  on  v  and  o,  we  note  that  its  deriv¬ 
ative  is 


do  EaT 

rear  _  ra  ■ — v  vo  -i 

dT  ~~  2(l-v)  ^l-o  (l-o)“ 


(62) 


Substituting  the  differential  equation  for  v,  we  find 


do  EoT  .  , .  ,  • 

rear  m  3v  r.  ( l-o)o-. 

dx  =  2( 1-v)  (l-o)?  L  3  J 


(63) 


Since  v  begins  (and  remains)  negative,  and  o  <  1,  the  derivative  of  the  stress 

remains  negative,  as  can  be  seen  from  Goodmans  computations  of  o(t)  and  the 

average  value  of  o,  which  may  be  bounded  above  by  the  average  recession  rate 

with  complete  removal,  o  <  4i/3. 
r  ave 

Thus,  within  che  limits  of  the  assumptions  made,  if  fracture  does  not 
occur  before  front  surface  melting  begins,  it  will  never  occur.  To  recap 
the  results  for  thin  slabs  (i  <  1)  in  terms  of  dimensionless  time 


kt 

T  =  U 


(64) 


and  intensity 


I  l 

i  =  -a- 

2kT  » 
m 

The  tensile  stress  on  the  rear  surface  is 


(65) 
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the  fracture  threshold  is  i  =  3oc,  and  fracture  times  for  two  hypothetical 

materials  are  as  shown  in  Figure  15. 

It  should  be  noted  that  all  of  these  results  are  for  the  boundary  con- 

3 

dition  considered  by  Laughlin  ,  i.e.,  no  bending  deformation  allowed,  and  no 
inplane  constraint  forces  allowed  to  develop.  Other  boundary  conditions  will 
produce  different  stresses,  as  will  be  considered  in  a  later  section. 
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Figure  IS.  Failure  Times  for  Thin  Slabs  Wleh  Melt  Retention. 
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2 


.  Thermal  Stress  in  Thick  Slabs,  with  Kelt  Retention. 

In  the  previous  section,  the  approximate  temperature  solution 
given  by  Goodman  and  Shea  was  used  to  predict  times  to  fracture 
in  thin  slabs  or  good  conductors  where  the  rear  face  i*  reached 
by  thermal  penetration,  and  hence,  the  maximum  stress  develops, 
before  front  face  melting  begins.  For  poor  conductors  at  high 
flux,  however,  the  front  force  will  melt  before  the  rear  surface 
heats.  In  such  a  case,  the  maximum  stress,  and  hence,  fractures, 
may  occur  after  melting  begins. 

We  therefore  require  a  solution  analogous  to  that  of  Goodman 
and  Shea,  applicable  to  the  thick  slab  geometry.  Let  us  assume 
the  one  dimensional  heating,  by  absorbed  intensity  la,  of  a  slab 
of  thickness  h.  We  denote  by  the  subscript  ra  the  thermal  proper¬ 
ties  of  the  molten  state,  and  the  depth  of  melting  by  s. 

On  0  <  x  <  s  3 1 Uffl  a  Um 


ons<x<s+5 

a2u  .  ajtj 

3x2  3 1 

The  boundary  conditions  are 


km  3x  "  la  x  -  0 


-  k 


3^m 

m  3x_ 


x«s 


+  kM 

+  K  3x 


X“  8 


T  d® 

pL  on  x 


Um"U“0onx"S 


3a 

s  3b 
4a 


k 


3U 

3x 


x«s+6 


0  on  x 


s+5 


U  ■  ■  Tm  on  x  -  s  +  6 

where  Tm  is  the  temperature  at  which  the  phase  transformation 

of  specific  heat  L  occurs,  and  temperatures  are  given  by 

T  ■  Tm  +  Ua  on  0  <  x  <  s 
T-Tn  +  U  ons<x<s  +  6 

The  initial  conditions  are  the  instance  of  front  face  melting. 


4b 

4c 


5a 

5b 


107 


U  ■  0  @  x  **  0  and  t 

a  ■  0  $  t  ■ 


U  -  -  TB  + 


(<5_  -  x) 

m 


t  ■  t, 


2k  - 


where 


6  *  6n  ■  /6k  ta  - 


2k  Tm 


«■  *V  UI  a  . 


and  the  initial  distribution  is  taken  as  the  approximate 
solution  prior  to  meeting. 

Following  Goodman,  we  choose  the  temperatures  distributions 
to  be  parabolic,  i.e., 

U  -  d(t)  +  xe (t)  +  x2f  (t)  .  ] 

Employing  the  boundary  conditions;  4a-c,  we  find 


U  ■  — ^  £s2  +  82o  -  2x  (s+o)  +  x2J 

or  T  -  T.  [(■-  '*  i  -  l]  -  Tm 

which  more  clearly  displays  the  physical  significance  of  the 
temperatures  ahead  of  the  melting  front  in  terms  of  the  two 
parameters,  the  melt  thickness  s,  and  the  penetration  depth,  6 
The  integral 


6  +  8 

0  -  /  Udx 
s 


'  T*  [  2  *] 


-f'. « 


will  prove  useful.  Physically,  it  is  a  representation  of  the 
additional  energy  required  to  bring  the  material  within  the 
penetration  depth  to  melting. 

He  also  choose  a  parabolic  temperature  in  the  molten  layer, 

Ua  n  a(t)  +  xb(t)  +  x2cCt)  14 

From  boundary  conditions  3a  and  4a, 

U*  ■  c(t)  Cx2-s2)  -  tx^s)  15 


108 


The  Integral 


6#  "  /  Undx  *  *'  *'  Y  ®3c(t)  16 

o  ® 

represents  the  energy  In  the  malt,  above  that  required  by  the 
phase  transformation.  At  t  «  0,  0O  «  0. 

We  now  require  the  satisfaction,  In  the  average,  of  the  heat 
conduction  equations,  (equations  1  and  2,}  through  averaging 
over  each  domain 


a 

/ 

0 

3U 

m  , 

j—  dx 

s  3U 

.  [  — ®  «JX 

>  3 1 

o 

17a 

and 

s+6  aiT 

/  K  37  dx 

s 

»+<S  an 

■  /  fr-* 

s 

17b 

Differentiating  13 

and  16,  we 

note : 

<36 

dt 

« 

«  (<5  +  8  J  U 

• 

|  -  a  U  + 

!  6+S  8 

s+d 

\  5T  dx 

18a 

and 

de 

m 

dt 

-  s  Um  + 

s 

s  3U 
atE  dx 

18b 

Hence 

3U 

*  _ 2. 

1C  1 

3  “  a 

de 

m 

1  Q_ 

*  3x 

s  sm 

dt 

and 

3U 

-  (C  r — 
3x 

d6  j.  r* 
dt  + 

+  «>  Tm 

19b 

1  8 

where  equation  4b  and  3a  have  been  applied.  Substituting  11  Into 
19b  leaves,  as  a  final  form  of  the  averaged  energy  equation  in  the 
solid,  that 


+2«  r-  ■  It  +  +  *>  T.  20* 

0  and  6  are  related  through  13,  thus 

e  -  -  2/3  Tm  6  21 
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Combining  IS  with  19a 


dO. 


2  •  c(t)  - 


22 


which  may  be  combined  with  16  to  give  the  a  final  averaged  energy 
equation  for  the  melt 


Ia«* 

2kT~ 


i  1  dl. 

-  e«  +  3  *2  7-  dF 

m 


23a 


We  may  also  eliminate  the  gradients  in  3b  to  obtain  the  equation 
for  the  energy  balance  on  the  melt  line  as 


I-  - 


dt 


2k  T. 


pLs 


24a 


Equations  20,  23  and  24  form  a  system  of  three  first  order,  non¬ 
linear;  coupled  equations  which  must  be  solved  in  order  to  find 
the  temperature.  Through  the  following  substitutions 

t  »  K  (t  -  tm)/i2  25 

M  ■  CpTm/L  26a 

v  ■  <b/k  26b 

i  -  Ia£/(2k  Tb)  26c 

we  arrive  at  the  same  dimensionless  variables  used  by  Goodman 


o  ■  s.  I 

f  -  0/(Tb  l) 

W  »  (tc  kn  9a)/(k  icm  T, 

The  differential  equations  become: 


£) 


27 

28 
29 


v  *  (1+y)  -  i  + 

w  "  3  (1  -  vw/o2) 

°  ■  v  { — r  +  T  7  '  i} 
0  i 


The  initial  conditions  Cx"  0)  are 

o  ■  0 
w  ■  0 

2 

V  "  ~  31 


20b 
2  3b 
24d 


30a 

30b 

30c 
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It  nay  be  demonstrated  that  if  t  1b  smell,  that 


o  *  6m1(1t)2 
w  ■  -  (6ui)2(iT)4 


V  .  _  J-  -  4iT  +  12wi(iT)2 


31a 

31b 

31c 


The  actual  temperature  distribution  in  the  solid  is  obtained  from 
12  as 


T  -  T  ( 

m 


32 


where  s  and  6  are  obtained  from  27,  13  and  2d.  Stresses  in  the 
solid  region  depend  on  the  mechanical  boundary  conditions  at  the 
edge  of  the  heated  region.  If  we  again  assume  that  bending  dis¬ 
placements  are  completely  constrained,  then 


E 

1 —v 


{Gi  -  oT> 


The  requirement  of  no  net  axial  force  leads  to 


Ea  1  f 
1-v  1  l-s  J 


s+6 


Tdx  -  T} 
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This  is  found  to  be 
EaT 


m 


{- 


1-v  3  (  i-s) 

EaT 


6  ,s  +  6  -  xv  2 


;  -  (: 


)2} 


m  /  6  .1 


for  x  <  s+6 


for  8+6  <  x  <  t 
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36 


”  1-v  'H-s'3 

The  maximum  tensile  stress  occurs  on  the  rear  face  (x  *  i)  and  is 

EaT 


rear 


2  (1-v)  l-o 


37 


valid  until  -  -j  v  <  l-o.  The  greatest  possible  stress  is 

EaT 
_ m 

°max  "  3 (1-v) 


For  small  time 


EaT 


rear 


(1 


fj£+  2it  -  6pi  (it) 2  +  2p (ir) 2  ) 


38 
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The  dimensionless  tine  at  which  salting  baglna  ia  (from  9,  25 
and  26c) 


and  tha  atraaa  at  that  time  ia 


EaT  , 

IB  1 

a  m  —  —  -  x- 

raar  1-v  3i 

Wa  nota  that  for  i  •  1,  tha  thraahold  of  applicability  of  thaae 
raaulta  for  tha  thick  alab,  that  tha  maximum  stress  predicted 
la  tha  aaae  as  for  the  thin  alab.  However,  we  also  note  an 
laportant  difference.  In  the  case  of  the  thin  slab,  we  found 
that  the  stress  reduced  after  the  onset  of  melting.  In  the 
thick  alab,  the  stress  will  continue  to  Increase  as  we  may  note 
from  Equation  38.  Hence,  a  numerical  solution  of  the  governing 
system  of  equations  is  necessary. 
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3.  Thermal  Stress  In  Melting  Thick  Slabs 
With  Instantaneous  Melt  Removal 

Let  us  now  consider  the  case  when  the  melt  is  removed  as  it  forms. 
Again,  two  cases  must  be  considered,  depending  on  whether  the  melting  be¬ 
gins  before  the  rear  surface  changes  temperature  (the  thick  slabs)  or  after 
(thin  slabs).  For  the  thick  slab,  we  assume  a  temperature  profile 


T  =  T  [(1  -  (^p)*] 
m  o 

for  s  <  x  <  s+6 

(1) 

T  =  0 

for  s+6  <  x 

(2) 

Here  s  is  the  instantaneous  location  of  the  front  surface  measured  from  its 
original  location.  This  temperature  profile  is  appropriate  only  after  front 
surface  melting  begins,  i.e«,  for 


2kT 


t  >  lm  =  6k 


_2h  * 


Again,  we  let 


8 


Udx 


2T  6 


m 


3 


(3) 


(4) 


where  U  =  T  -  T^.  Requiring  the  satisfaction  of  the  differential  equation 
averaged  over  the  instantaneous  thickness  leads  to 

- "  s  i  ■  ecp  r  if  "x  ■  •  *  Vs + <5> 

1  s  r  s 

The  melting  line  coincides  with  the  front  surface  and  an  energy  balance  leads 
to 
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ax  i  • 

I  +  k  —  =  pLs 

a  3x 

s 


(6) 


The  thermal  gradient  on  x  =  s  may  be  evaluated  from 


3T 

3x 


2kT 


m 


(7) 


Application  of  Equation  ( 4  )  and  the  elimination  of  the  gradients  leads  to 
two  equations 


•  T  *  V 

°Ls  =  Xa  +  T  -T 


(8) 


PC  -r  —  p(L  +  C  T  )s  =  -I 
p  i.  pm  a 


(9) 


In  terms  of  the  same  dimensionless  parameters  introduced  previously,  we  find 


dv  .  8  ..  >1 

—  =  +  -r  (1  +  u)  - 

dr  3  v 


(  J) 


£  ■  *  £ 


(11) 


Again,  the  initial  conditions  (t  -  t  =0)  are 

IQ 


o  =  0 


(12a) 


v  =  - 


31 


(12b) 


In  this  case,  we  note  the  first  equation  to  be  uncoupled.  The  substitution 


,7  =  C[1  +  g(t)] 


(13) 
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is  found  to  be  convenient.  With 


C  = 


3  ui 

2  (1+w) 


(14) 


and  the  change  of  independent  variable 


y  =  f  (1+u)  c*t  =  yf-  M«i* r* 


(15) 


where 


T 


I 


=  T  -  T 

in 


(16) 


then 


dy 


=  -  g(l  +  g)* 


with  initial  condition 


8(0)  =  J 


This  equation  may  be  integrated  to  yield 


y  =  -  -p-  +  In  (-A±2_) 

1+u  1+g  g(l+u) 


Substitution  into  Equation  (11)  then  yields 


(17) 


(18) 


(19) 


3pio  =  u 


)  +  In  ( 


_i±B_ 

g(l+u) 


) 


(20) 


This,  together  with 


vi 


-2(1+11) 

3m( 1+g) 


(21) 
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enables  the  solution  to  be  computed  in  terms  of  the  parameter  g,  which  is 
in  turn  related  to  y,  and  hence  t  through  (19)  and  (15).  Pertinent  values 
of  solution  parameters  are  presented  as  Table  I. 

For  the  same  boundary  condition  as  was  considered  in  the  previous  sec¬ 
tion,  i.e.,  no  bending  displacement  or  inplane  force,  we  find 


a 


ee 


Tdx  -  T] 


(22) 


Of  particular  interest  is  the  maximum  tensile  stress,  which  again  develops 
on  the  rear  face,  x  =  Z,  where  T  =  0  for  times  over  which  this  thick  slab 
solution  appHes.  Hence, 


EaT 

_  _ m 

°niax  ~  2(l-v) 


(23) 


In  this  case,  a  power  series  solution,  valid  for  small  time,  shows  the  stress 
at  the  rear  face  to  increase  after  front  face  melting  begins,  since  both  v 
and  a  increase  in  magnitude. 

Values  of  a  dimensionless  stress 


a 

max 


o  (1-v) 
max _ 

EaT 

m 


-  -  (— ) 
2  U-o' 


(24) 


may  be  computed  for  any  i  >  1. 

It  is  to  be  recalled  that  the  temperature  distribution  becomes  invalid 
when  the  tht-mal  penetration  front  reaches  the  rear  face,  or 

a+2=l  (25) 


This  condition  occurs,  for  any  i,  at  a  time  y  when  i  =  i£,  as  determined  from 
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TABLE  I.  SOLUTION  PARAMETERS  COMPUTED  FOR  \i=2. 


& 

1 

-Vi 

io 

i  =io-riv 
c  2 

0 

c 

5 

0 

.66667 

0 

1.0 

45 

.04847 

.68965 

.000415 

1.0349 

.000401 

4 

.10653 

.71430 

.0018817 

1.0733 

.00175 

35 

.17724 

.74074 

.004848 

1  >  1 160 

.00434 

3 

.26516 

.76923 

.010005 

1.1639 

.00860 

25 

.37749 

.8000 

.018472 

1.2185 

.01516 

2 

.52648 

.83333 

.032192 

1.2822 

.0251 

15 

.73537 

.86955 

.05493 

1.3593 

.0404 

10 

1.0569 

.90909 

.095335 

1.4590 

.065 

07 

1.3604 

.93460 

.137427 

1.5394 

.0893 

05 

1.6602 

.95238 

.18146 

1.6100 

.113 

03 

2.1333 

.97087 

.25415 

1.7105 

.149 

02 

2.5195 

.98039 

.31534 

1.7859 

.177 

01 

3.1931 

.9901 

.42437 

1.9095 

.222 

001 

5.4768 

.999 

.802182 

2.302 

.348 

0001 

7.7786 

.9999 

1.18535 

2.685 

.441 

00001 

10.081 

.99999 

1.56906 

3.0690 

.511 

.02  x  10  7 

15.667 

1.0 

2.5 

4.0 

.625 

k.9  x  10~10 

20 

1.0 

3.2222 

4.7222 

.682 

:  .53  X  lo"10 

21.6667 

1.0 

3.5 

5.0 

.70 

:  .24  x  lo"14 

30 

1.0 

4.8889 

6.3889 

.765 

.85  x  10-18 

39.6667 

1.0 

6.5 

8 

.8125 

.014  x  10-18 

40 

1.0 

6.5556 

8.0556 

.813 

:.7  X  10"23 

51.667 

1 .0 

8.5 

10 

.85 
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the  appropriate  column  in  Table  I.  At  that  instant,  the  stress  is  a  maximum 
and  has  value 


o 

max 


1 

3 


(26) 


independent  of  i."  Thus,  the  maximum  stress  obtainable  is  independent  of  the 
dimensionless  intensity,  and  u  as  well,  but  the  time  required  to  generate 
the  stress  will  be  seen  to  depend  on  the  intensity.  The  column  labled  a  _ 
in  Table  I  denotes  the  fraction  of  the  slab  which  has  melted  at  the  instant 
when  the  stress  achieves  the  maximum  value. 

Figure  16  presents  the  time  history  of  the  evolution  of  the  rear  face 
tensile  stress,  for  various  values  of  the  dimensionless  flux.  In  Figure  17. 
the  failure  time  as  a  function  of  intensity  is  given  for  materials  of  various 
ultimate  strength,  o^,  where 


A 

o 

u 


Vtt(1-v) 

EaT 

m 


(27) 


One  notes  what  at  first  would  appear  to  be  an  anomolous  increase  in  fracture 
time  with  intensity.  This  arises  because  higher  intensities  confine  the 
penetration  depth,  thus  less  thermal  force  develops.  Further  increases  in 
intensity  cause  more  rapid  melt  thru,  thus  reducing  the  load  bearing  area. 
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Figure  16.  Stress  Time  History  on  Rear  Surface  of  Thick  Slab  With 
Melt  Rcraova 1 . 
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Figure  17.  Failure  Times  for  Thick  Slabs  With  Melt  Removal 


V.  Summary  and  Discussion 


Two  mechanisms  have  been  identified  as  significant  in  the  generation  of 
stress  in  solid  elastic  materials  subjected  to  high  intensity  laser  radiation. 
The  absorption  of  some  fraction  of  the  incident  energy  was  found  to  produce 
thermal  stresses  of  significant  magnitude  ,  but  these  stresses  are  character¬ 
istically  compressive  in  the  region  heated  by  the  beam.  In  a  ductile  material, 
such  stresses  may  produce  failure  in  a  shear  mode,  but  this  is  less  likely  to 
occur  in  a  brittle  material.  Rather,  it  is  the  significant  tensile  stresses 
found  outside  of  the  heated  region  which  are  believed  to  be  the  more  significant 
in  the  generation  of  failure.  Radial  conduction  of  heat  was  found  to  diminish, 
but  not  remove,  these  tensile  stresses. 

Such  thermal  stresses  are  generated  by  both  CW  and  pulsed  lasers.  In  the 
latter  case,  the  duration  of  the  pulse  and  the  time  between  pulses  is  usually 
small  compared  to  the  thermal  diffusion  times;  thus  a  repeated  pulse  may  be 
treated  as  an  energy  source  of  the  same  time  averaged  intensity.  This  approx¬ 
imation  is  especially  appropriate  in  the  case  of  poor  conductors,  such  as 
glasses  and  ceramics. 

In  addition  to  the  stresses  developed  by  the  absorption  of  energy  by  the 
solid,  the  absorption  of  energy  in  the  medium  adjacent  to  the  solid  produces 
an  intense  pressure  field  through  the  mechanism  of  LSD  or  LSC  waves.  The 
mechanical  response  of  an  elastic  plate  to  such  loadings  was  also  considered. 
This  mechanism  is  inherently  of  much  greater  significance  in  the  case  of  pulsed 
lasers  because  of  the  higher  peak  intensities.  Consideration  of  the  response 
to  a  pulse  showed  that  the  response  of  a  short  pulse,  acting  over  a  small 
portion  of  the  plate,  can  be  accurately  determined  only  if  a  large  number  of 
the  plate  modes  are  considered  in  the  response.  While  the  peak  stresses  occur 
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at  the  center  of  the  spot,  the  stresses  elsewhere  In  the  plate  remain  signifi¬ 
cant. 

The  stresses  produced  by  these  two  mechanisms  are  such  that  superposition 
occurs  in  the  case  of  a  repeated  pulse  laser  of  high  Intensity.  The  thermal 
stresses  due  to  the  time-averaged  intensity  increase  with  time.  At  later  times, 
the  peak  values  are  to  be  found  at  greater  distances  from  the  heated  area.  The 
oscillating  stresses  resulting  from  the  front  surface  pressures  will  diminish 
in  magnitude  with  distance  from  the  spot  center.  While  oscillating  at  the  "rep" 
frequenc  '  of  the  laser,  the  amplitude  of  succeeding  peaks  will  be  relatively 
constant.  Thus,  the  superposition  of  the  two  stress  generating  mechanisms  may 
well  produce  failure  after  a  number  of  pulses  have  occurred.  Failures  occurring 
in  the  first  few  pulses  may  be  expected  to  be  due  primarily  to  the  oscillating 
pressure  fieM  and  therefore  at  or  near  the  spot  center  (barring  stress  raisers 
at  the  supports).  Failures  which  occur  after  a  number  of  pulses,  on  the  other 
hand,  r  y  1  i  expected  to  be  due  to  the  development  of  the  tensile  thermal  stress 
outside  o_  che  heated  region.  Thus,  the  failure  may  originate  outside  of  the 
heated  region. 

Other  fai  ure  mechanisms  were  given  some  consideration.  Tensile  stresses 
were  found  ;•*  be  present  in  the  interior  at  early  times,  but  the  reported 
significance  of  surface  flaws  as  failure  origins  is  not  accounted  for  by  these 
stresses.  Partial  melting  under  the  spot  may  be  of  significance.  A  preliminary 
investigation  of  this  possibility  was  conducted,  with  the  scope  limited  by  the 
mechanical  boundary  condition  considered  at  the  edge  of  the  melting  region. 
Briefly,  the  failure  mechanism  deserving  some  further  exploration  is  the  fol¬ 
lowing.  As  melting  occurs,  the  central  portion  of  the  plate  becomes  thinner. 

The  compressive  forces  generated  by  the  surrounding  unmelted  portion  of  the 
plate  will  then  have  a  resultant  nearer  to  the  melted  surface  than  to  the  rear, 
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thus  generating  the  proper  moment  to  produce  rear  surface  tension.  While  such 
e  mechanism  is  conceivable,  it  is  not  likely  to  be  of  significance  in  a  ceramic 
because  of  the  long  times  required  to  produce  the  necessary  melting  (or  thermal 
softening)  necessary  for  the  process  to  occur. 

While  the  possibility  of  crack  (flaw)  growth  under  the  oscillating  stress 
produced  by  the  fluctuating  front  surface  pressure  cannot  be  precluded  as  a 
producer  of  the  observed  failures  after  a  number  of  pulses,  thermal  stresses 
of  increasing  amplitudes,  upon  which  the  fluctuating  stresses  are  superposed, 
would  seem  to  be  a  more  probable  cause.  A  carefully  designed  and  properly 
instrumented  series  of  experiments  would  establish  the  appropriateness  of  this 
conjecture. 
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